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ABSTRACT 

Ground motion characteristics of the Himalayan earthquakes have 
been studied from the 198 horizontal component records and 99 vertical 
component records from six earthquadces in the region that are now 
available. The characteristics studied are the peak parameters, duration 
of strong shaking, shape of elastic response spectrum and ductility 
reduction factor. It is found that mean A/V ratio is quite high in the 
region and does not vary much with the distance from the shock-source. A 
study of shape of elastic acceleration respcaise spectra revealed that: 
(1) shape varies somewhat insignificantly with the distance, (2) shape 
of spectra in the intermediate period range highly depends on earthquake 
magnitude, (3) there is a urgent need to revise the shape of response 
spectra recommended in 13:1893-1984, and (4) shape of "standard 
spectra" suggested by AERB for rock sites is comparable with the 
raean-plus-one-standard-deviation shape from the data. Both, the ratios 
of peak parameters and elastic response spectra indicate that with 
distance the attenuation of waves of different frequencies is not so 
. different as elsewhere. Based on the data, shape of ductility reduction 
factor spectrum to be used in design is also proposed. 
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1.1 GENB»L 

Earthquake-resistant design of structures requires estimate of the 
ejqpected ground motions at the site in future earthquakes. In order to 
develqp ability to reasonably assess the expected ground motion due to 
future earthquakes, it is essential to study the characteristics of 
strong ground motions recorded in the past earthquakes. 

Design of engineering structures is usually governed by building 
codes and standards of design practices. Hie recomnnendations in these 
codes are often based on the analysis of records from major earthquakes. 
For example, the seismic design spectrum, specified in the codes, takes 
into consideration characteristics of recorded ground motions. As more 
records become available, these recommendations are updated and improved 
for better earthquake-resistant design. Unfortunately, until recent 
years, time history records of strong ground motion in the Indian 
earthquakes were almost non-existing. As a result, for instance, the 
shape of design spectrum in the Indian seismic code (IS: 1893 - 1984) is 
based mainly on some earlier earthquakes in California. 

Since 1986, a number of strong motion records have become available 
from earthquakes in Himalayas. These need to be studied and analysed to 
understand the peculiarities of ground motion in the Himalayan region. 


1.2 SCOPE OF THE STODT 

In the present work, strong motion instrument data from Himalayas, 
recorded during the period 1986-1991, are examined. In all, 198 
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horizontal carponent and 99 vertical ccmponent records from six 
different earthijuakes with magnitude varying from 5.2 to 7.2 and 
epi central distance in the range of 5 knt to 320 km are now available. 
These have been analysed to understand the characteristics of 
earthquake ground motion in the Himalayas. The characteristics studied 
are the peak parameters, duration of strong shaking, shape of respcmse 
spectrum, and ductility reduction factor. Results obtained are compared 
with those from published studies based on earthquakes elsei^ere and 
with currently adopted design practices. 

1.3 ORGANIZATION 

Apart from this introductory chapter, the thesis comprises of five 
more chapters to cx3ver the present work. Chapter II gives details about 
the strong motion records used in this study. Chapter III presents the 
characteristics of the data in tenns of ratios of peak ground motions 
and duration of strong shaking. Chapter IV presents the study on shape 
of response spectra obtained from the data. Chapter V presents the sttdy 
on ductility reduction factor for elasto-plastic model. Summary and 
conclusions from the presaat work are given in Chapter VI. Appendix A 
gives the details of the records used in the study with regard to 
name of stations triggered in each earthquake, epicentral distances, 
peak ground acceleraticn, peak ground velocity, peak ground 
displacement, duration of strong shaking (Trifunac and Brady 1975), and 
spectrum intensity (Housner 1959) for 5% damping. 



CHAPTER II 

STRONG MOTION RECORDS 


2.1 raraODOCTIOK 

Ground motion during an earthquake is measured by a seismograph 
which is essentially a single degree of freedom system with its mass, 
stiffness, and damping characteristics appropriately designed. 
Instruments with large natural period (flexible oscillator) usually 
record ground displacement; these being flexible are very sensitive and 
can record ev^ the slightest motions in ground. Such instnxnents, 
notmally termed as seismographs, are useful for recording very small 
earthquakes or for recording large earthquakes taking place at large 
distance away from the instrutnent. However, these go off-scale in case 
of a strcxig shaking, which is of primary interest for engineering 
a^pplications. On the other hand, instruments with a stiff oscillator 
(large natural frequency) measure ground acceleration and are not 
sensitive to weak ground motiais. Such instruments, often tenned as 
strong motion accelerographs, provide the realistic levels of ground 
acceleration that an engineering structvire nay be siA)jected to during a 
strong shaking. In normal operations, these instnxnents lie in a passive 
state, and get triggered, to operate, when the grovmd motion exceeds a 
threshold value set for the instnanait. 

The acceleration versus time records obtzdned from such instnanents 
are called accelerograms. Usually, an accelerogram is in the form of an 
analog record, such as a photographic film. However, in recent years, 
digital instrunents have also become available. The analog records 
obtained from the accelerograph films have to be digitized at closely 
spaced time intervals. The digital time history, obtained directly or 
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from digitizing the analog record, is then corrected through appropriate 
filtering to control the digitization errors and baseline distortion; 
also certain accelerograph transducer corrections are applied. The 
records thus processed are termed as the "corrected" accelerograms. The 
corrected accelerogram can then be integrated with respect to time to 
obtain a velocity versus time record, which can be frurther integrated to 
obtain a displacement versus time record. 

2.2 CaaRACTHlISKriCN OF SraCNO GROOND MOTION 

Pig. 2.1 shows some typical accelerograms obtained from epicentral 
region of different earthquakes. A visual comparison of these records 
suggests that the characteristics of earthquake ground motion may vary 
significantly with regard to the peak acceleration, number of 
significant peaks in the time record, and duration of shaking. 

Pig. 2. 2 shows the acceleration, velocity, and displacement time 
histories derived from a typical accelerogram. Not surprisingly, the 
integration of record changes its frequency content; ^diile the 
acceleration record has large high frequency content, the displacement 
record shows motion of high period, and the velocity record shows an 
intermediate frequency content. A structure is expected to achieve 
resonance (or multifold increase in response) when frequency of ground 
motion matches with natural frequency of the structure. Hence, 
structiures with different natural period are expected to respcnd 
differently in the same earthquake. 

The earthquake ground motion, which is somewhat randan with time, 
can be characterized by several parameters depending on the application. 
•Hie parameters of most interest are the peak parameters (peak ground 
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acceleraticm, velocity, or displacanent) , duration of strong shaking, 
and its frequency contoit. The frequency content is somewhat more 
difficult to quantify and is usually represented through a Fourier 
spectnm or a response spectrum. 

Traditionally, in the earthquake- resistant design, peak ground 
acceleration has been the single most frequently used parameter for 
describing the severity of ground motion. A rigid structure subjected to 
a given earthquake ground motion experiences a total earthquake force 
equal to its mass times the peak ground acceleration. Ibus, peak ground 
acceleration gives a good physical feel of the ground motion and is a 
ccaivenient rule of thimb for scaling different accelerograms. 1510 
difficulty with peak ground acceleration is that, sometimes, it gives a 
very misleading indication of the earthquake ground motion. For 
instance, an earthquake of rather low magnitude (say M = 4 or 5) may 
give a rather large peak ground acceleration, even though the overall 
motion may not be very damaging to structures. 

The peak ground velocity is associated with waves of intermediate 
frequQQcies (mostly betweoa 0.3 to 3.0 Hz), which is also the range in 
which most civil engineering structures lie. Hence, peak ground velocity 
is a much better indicator of the damage potential of a ground motion 
than the peak ground acceleration. Peak ground displacement is 
associated with waves of low frequencies (0.1 to 0.3 Hz); in this 
frequency range very few structures lie, hence this parameter is not 
very useful for engineering applications. However, peak ground 
displacement is an important parameter for seismological studies, 
because, amplitudes of ground displacement are propagated more 
coherently than amplitiides of ground velocity and acceleration. This due 
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to the fact that their low-frequency spectral carposition is not very 
sensitive to scattering by small geologic inhonnogeneities. Mso, peak 
ground displacement gives a scmewhat better indication of the fault 
displacement . 

2.3 STRONG MOTION mSISSCIlENTATION IN 

The arc-shaped Hinalayan belt extending from Pamirs in erstvdiile 
USSR to Arkans in Burma has long bea^ regarded as a highly seismic area 
in the world. The Himalayan belt is believed to have formed due to 
collision of Indo-Australian and Eurasian plates in past several mil lion 
years. The area has experienced several destructive earthquakes of 
magnitude greater than 8.0. The notable among these earthquakes are the 
Assam earthquake of 1897 (M = 8.7)/ the Kangra earthquake of 1905 (M = 
8.6), the Bihar-Nepal earthquake of 1934 (M = 8.4), and the IVssam-Tibet 
earthquake of 1950 (M = 8.7) (Fig 2.3). Munerous eairtbpakes of 
magnitudes ranging between 6.0 and 7.0 have also occurred in the region. 

An intematicnal workshop on strong motion earthquake instnment 
arrays held in Honolulu, Hawaii, in May 1978 (Iwan 1978) identified the 
area around Shillong in the north-eastern India as the highest priority 
site for deployment of a strcxig motion array. Subsequently, three strong 
motion arrays, listed below, have been d^loyed in the Himalayan region: 

(a) An array of fifty strong motion accelerographs (SMA-1 of m/s 
Kinemetrics, USA) is located in north-west Himalayas in the state of 
Himachal Pradesh. This array is identified as Kangra array. 

(b) An array of fortyfive SMA-1 is located in north-east Himalayas in 
the states of Assam and Meghalaya. This array is identified as 


Shillong array. 
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(c) An array comprising of forty SMA-1 instruments is located in the 
central Himalayas in the state of Uttar Pradesh. This is identified 
as Uttar Pradesh array. 

Fig. 2.4 shows the location of these arrays. There are also plans 
to instrument the Himalayas of Bihar-Nepal border region. Thus, the 
Himalayan zone is now better instnmented to record the strong-motion 
data from future earthquakes. 

Besides the above three arrays, the region also has a number of 
indigotious strong motion accelerographs and structural response 
recorders (BQ 88-16 1988, EQ 93-16 1993); these have also provided 
records of strong ground motion. However, the data from these 
instninents has not been considered in the present study. 


2.4 STOGNG MOTION DATA USED IN TOE PRESENT STOD? 

The Kangra array registered an event of magnitude 5.5 on l^ril 26, 
1986. This is the only significant event recorded so far by this array. 
TIm aiillong array has registered four events during 1986-1988, while 
the Uttar Pradesh array has registered one event in October 1991. The 
details of these events with regard to their date, magnitude, maximum 
peak ground acceleration, and the range of epicentral distance are 
presented in Table 2.1. This table also gives the abbreviation used in 
rest of this study for describing different events (e.g., "Sept” for the 
earthquake of September 10, 1986). The magnitudes indicated in table are 
as per estimates of the United States Geological Survey (USGS); these 
differ somewhat from those by the India Meteorological Department (IMD). 
Thus, the present study is based on 99 sets of strong motion records 
(198 horizontal components and 99 vertical components) obtained from six 
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earthquakes (magnitude ranging frcxn 5.2 to 7.2) that have ooourred in 
three different regions of the Himalayas. The epicentral distance for 
these records vauries from about 5 km to 320 km. Detailed description on 
site geology of stations is not available. However, Chairirashddbaran and 
Das (1990) indicate that the stations are located on "reasonadily firm 
ground", "soil", or on "soft sedimentary rock". In the absence of more 
specific information about the local soil conditions, all stations are 
assumed to have almost similar groxmd conditions, i.e., the recscards have 
not been classified on the basis of soil conditions. Moreav&c, vdiile 
comparing with other data, the stations are assuned to be at rock sites. 

Figures 2.5 to 2.10 show the location of instrxinents triggered 
during the different earthquake events. Also shc»jn in these figures is 
the location of epicenter. Figs. 2.5 to 2.9 show three locations of 
epicentre: El is the estimated location by the IMD, E2 is that by the 
uses, and E3 is the location st^gested by Chandrashekharan and Das 
(1990) based on records from the array. For "April" (Fig. 2.5), revised 
estimate of epicentre by IMD, E4, is also shown. For each of the events. 
Appendix A gives the name of the stations triggered and three different 
estimates of epicentral distance for each station. Also presented in 
Appendix is the peak ground acceleration, peak ground velocity, peak 
ground displacement, spectrum intensity (Housner 1959) for 5% damping 
and duration of strong shaking (Trifunac and Brady 1975) for each of the 
records used in the study. The different estimates of epicentre in these 
earthquakes vary by as much as 68 km for "Sept" and 85 km for "Feb". 
Such large variation niles-out the possibility of using epicantral 
distance as a criteria for classifying the records for low, medium, and 
large epicentral distances. Moreover, for most of these events, neither 
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an indicaticjn of n 5 )txire zone (based on aftershock data) nor the 
isoseisnials are available. Thtis, the distance of recording station from 
the zone of energy release also cannot be used as a criteria for 
classifying the records. However, in this stxidy, epicentral distance 
refers to the ^icentral location sijggested by Chanshekhran aM Das 
(1990) based on the analysis of strong motion records. The location of 
^icentre for the "Oct” event (Fig. 2.10) is fairly accurate. 

To account for the variation in epicentral distance, each set of 
records has been classified into one of the three groups on the basis of 
the value of larger horizontal peak ground acceleration (PGR). Records 
with larger horizontal peak ground acceleration less than 0.05g are 
classified into the low PGR group, whereas those with greater than O.lg 
PGR are classified into the high PGR group; records with larger 
horizoital peak ground acceleration between these two limits are 
classified into the medium PGR group. Thus, once the larger horizontal 
PGR from a set falls in one gro^p, all the three components of that set 
are included in that groifl?. The limits of PGR for classification were 
chosen so that in each group a reasonable number of accelerogranns could 
be included. Table 2.2 shows the surmary of ntnher of records included 
in each of the groups from different events. 
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Table 2.1 Summary of the earthquake events ccxisidered in IMs study. 


Array 

region 

Date of 
event 

Abbreviated 

as 

M 

NS 

R (km) 

PGA (g) 

Min. 

Max. 

Horz 

Vert 

Kangra 

April 26, 
1986 

April 

5.5 

9 

5 

29 

0.250 

0.082 

Shillong 

S^teniber 10, 
1986 

Sept 

5.2 

12 

15 

181 

0.138 

0.062 

May 14, 

1987 

May 

HR 


46 

179 

0.086 

0.005 

February 6, 
1988 

Feb 

5.8 

18 

35 

176 

0.114 

0.103 

August 6, 

1988 

Aug 

7.2 

33 

100 

324 

0.343 

0.180 

Uttar 

Pradesh 

October 20, 
1991 

Oct 

6.6 

13 

25 

151 

0.310 

0.294 


M : Magnitxide of earthquake 
NS : Number of stations triggered 
R : E);>icentral distance (approximate) 
PGA : Peak ground acceleration 
Horz : Horizontal component 
Vert : Vertical component 


Table 2.2 Sunnary of nuriber of records in the three PGA groiqps. 


Array 

region 

Event 

High PGA 
group 

Medium PGA 
group 

Low PGA 
group 

Horz 

Vert 

Horz 

Vert 

Horz 

Vert 

Kangra 

April 1 

10 

5 

4 

2 

4 

2 

Shillong 

Sept 

6 

3 

8 

4 

10 

5 

my 

0 

0 

12 

6 

16 

8 

Feb 

2 

1 

12 

6 

22 

11 

Aust 

28 

14 

32 

16 

6 

3 

Uttar 

Pradesh 

Oct 

8 

6 

12 

4 

6 

3 
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CHAPTER III 


RATIOS (F PEAK GROUND MOTION AND 
DURATION OF STRONG SHAKING 


3.1 nmanociim 

Given an ^timate of horizontal peak ground acceleration (PGA) at a 

given site, the values of peak ground velocity (PGV) and displacement 

(PGD) are often obtained by using average values of ratio of peak 

2 

parameters. Two ratios, A/V and AD/V (A: peak ground acceleration, V: 
peak ground velocity, D: peak ground displacement) are most often used 
for this purpose. Similarly, vertical peak ground acceleration is 
usually estimated by using an average value of the ratio A^/Ag (A^ : 
vertical peak ground acceleraticai, Ag : horizontal peak ground 
acceleration). A number of studies are available idiich indicate the 
average value of these ratios for a given type of site condition, based 
on ground motions recorded elsewhere in the world (e.g. , Mohraz et al. 
1972, Mdiraz 1976, Mohraz 1991). 

Seismological studies show that attenuation of seismic waves is 
frequency dependent; high frequency waves attenuate more rapidly with 
distance than moderate or low frequency waves. Since, peak ground 
acceleration, peak ground velocity, and p^k ground displacement are 
associated with waves of different frequaicies, these ratios of peak 
parameters happen to be good indicators of ground motion characteristics 
at different sites. For example, ground motions experienced near an 
earthquake source are expected to have higher value of A/V than ground 
motions at a large distance from the source of seismic energy release. 
This is so because the high frequency waves (and hence the PGA) 
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attenuate faster with distance than the intermediate frequency waves 
(and hence the PGV) . 

2 

It can be seen that for harmonic oscillations, AD/V = 1.0. With 
distance, the ground motion tends to become harmonic with 

ever-increasing period; hence as epi central distance increases, the 

2 

AD/V ratio must reduce and approach to unity (Newmark and Rosenblueth 
1971) . Also, certain response spectirum characteristics can be correlated 
to the AD/V^ ratio. 

This chapter presents a study of the peak parameter ratios and 
duration of strong shaking from strong motion data of the six Himalayan 
earthquakes . 

3.2 REVIEW OF TOE PAST STODIES 
3.2.1 A/V Ratio 

Mohraz et al. (1972) conducted the first study on V/A ratio. Tliey 
obtedned the average values of V/A from a data set, mostly consisting of 
records obtained from the San Fernando earthquake of 1971. They 
recotmvended the value of V/A ratio as 48 in/sec/g (A/V ratio as 8.1 
sec ^) for firm ground, and 32 in/sec/g to 36 in/sec/g (A/V ratio in the 
range 11.7 sec ^ to 12.1 sec ^) for rock sites. 

Mohraz (1976) studied the variation of V/A ratio with different 
site conditions using 54 earthquake records obtained in 16 seismic 
events in Califomnia. The records were classified into four categories 
depending on the geological conditions at site. For each site category, 
the ratios were obtained for three groups; the first group consisted of 
horizontal coirponent records with larger PGA, the second of horizontal 
corponent records with smaller PGA, and the third of vertical component 
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records. Table 3.1 sxitmarises these results. It can be seen that A/V 
ratio is higher for ro<dc sites than for alluviim. Also, A/V ratio for 
horizontal conpcxiait groi?) with larger PGA is higher than that for the 
lower PGA grou®*. A/V values for the vertical ccrnponent are close to 
those for horizontal cornponent grov^ with larger PGA. 

Variation of V/A ratio (and hence A/V ratio) with distance from the 
njptiared surface (R) has been studied by Mohraz (1991) from records of 
the Lam Prieta earthquake of 1989 (Table 3.2). It is seen that at rock 
sites A/V ratio for near-field records (R < 20 km) can be more than 
twice that for mid-field records (50 km > R > 20 km). At alluvium sites, 
the difference is about 30 %. However, as the site distance increases 
beyond 50 km, A/V ratio shows a slight increase. 

Using 36 records from rock sites, Zhu et al. (1988) studied the 
effect of ground motioaos with low, normal, and high A/V ratio on the 
damage of single degree of freedom stiffness-degrading systems. They 
demonstrated that if a system is subjected to groiand motions with same 
PGA as that of the design spectrum, the ground motions with low A/V 
ratio can cause much more significant peak inelastic deformation, 
stiffness deterioration, and hysteretic energy dissipation than is the 
case for ground motions with high A/V ratio. This is consistent with the 
general ejqpectation that given two ground motions with the same PGA, the 
motion with higher PGV will be more damaging to the structiires. 

Tso et al. (1*992) have examined the engineering implications of A/V 
ratio on the ground motion characteristics. A data set consisting of 45 
records obtained from rock sites during 23 events having nagnitrde 
ranging between 5.2 to 8.1 was used. The records were classified as low, 
normal, or high A/V when this ratio was less than 0.8 g/m/s ( = 7.9 
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sec ^)/ between 0.8 g/m/s and 1.2 g/m/s ( = 11.8 sec and more than 
1.2 g/m/s. The study concluded the following: 

- Ground motions in the vicinity of small or moderate earthquake usually 
have high R/V ratio, whereas those distant from large earthquakes have 
low R/V (Fig. 3.1). 

- Ground motions with high R/V ratio xjsually have seismic energy in the 
high frequency range, whereas those with low R/V ratio have energy in 
low frequency range. It was seen that ordinates in acceleration 
spectra for periods longer than about 0.5 sec for low R/V records 
could be more than twice as high as those for high R/V groiq? (Pig. 
3.2). 

- Peak ground velocity correlates well with the spectnm intensity 
(Housner 1959) over earthquake records having drastically different 
R/V ratio. 

3.2.2 M)/V^ Ratio 

The RD/V^ ratio was first introduced by Esteva (1969) . The ratio is 

primarily used to estimate PS) from known PG3\, and PGV. He suggested an 

2 

empirical equation for estimating AD/V for a known epicentral distance 
(R in km): 

M 1 4 . 400 

rO.6 

2 

Newnark and Rosenblueth (1971) suggest that choice of AD/V to evaluate 
P® is preferable as compared to D/A or D/V ratios beca\;ise the 
integrated velocities and displacements are sensitive to baseline 
adjustment and large errors may accumulate during the integration. By 
selecting AD/V , the error which may arise is distributed among the 
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ground motion parameters. According to Newmark and Rosenblueth, for most 
earthquakes of practical interest AD/V^ ranges ' from approximately 5 to 
15. 

McAoraz et al. (1972) reconmend an "adequately conservative" value 
of 6.0 for AD/V^ ratio. Further, they observed that AD/V^ is a measure 

of sharpness or flatness of pseudo velocity spectnxn, i.e., higher the 

2 

value of AD/V ratio, flatter will be the pseudo velocity spectrum. 
Table 3.1 clearly shows that AD/V^ is higher for rock sites than for the 
alluvion; however, this is not so obvious from the Lena Prieta 
earthquake data (Table 3.2). Also, Tetble 3.2 does not show any specific 
trend in the valtie of AD/V^ ratio with site distance from the ruptured 
surface. 

3.2.3 Ratio 

Table 3.3 shows the mean and 84.1 percentile values of the vertical 
to larger horizontal PGA ratio (Ay/A^) edatained by Mohraz (1976). For 
rock sites the mean value is around 0.52 while the 84.1 percentile value 
is around 0.69. He dDserves that "The ratio of the vertical to larger 
peak horizontal acceleraticai of 2/3, which is curraitly employed in 
seismic design corresponds to 84.1 percentile ratio and not the mean. 
Although this value is conservative, its use has been justified to 
account for variations greater than the median and mean". The Indian 
code (15:1893-1984) reccranends that the ratio of vertical to horizontal 
seismic coefficient be taken as 0.5. 

3.3 RESULTS ERCM THE DATA 

For each of the six events, mean and mean-plus-standard deviaticai 
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values of the h/V and AD/V^ ratios for high, medium, and low PGA groi:i»s 
have been computed; these are siamarized in Table 3.4. The ratio for 
vertical PGA to the larger horizontal PGA (\/^) for different events 
is given in Table 3.4. The trends for the Himalayan earthquakes from 
these tables are discxassed here. 

3.3.1 A/V Ratio 

Table 3.4 clearly shows that the A/V ratio for the Indian 

earthquakes is unusually higix. For the horizontal compmeats, the mean 

-1 -1 

value for different events ranges from 13.3 sec to 32.5 sec with an 
overall mean value of 20.2 sec These values are to be ccnpared with 
mean values of 12 sec ^ to 15 sec ^ for rock sites suggested by 
different studies of Mohraz based on data for elsewhere. For instance, 
Mohraz (1972) suggests the value as around 12 sec"^; Table 3.1 indicates 
a value of 14.5 sec and the data of the Lama Prieta earthquake (Table 
3.2) indicates a mean value of 14.8 sec ^ for near-field (and 
significantly lower values for site distance exceeding 20 km). 
Similarly, the data sttxJied by Tso et al. (1992) pertains to epicentral 
distance ranging from 4 km to 379 km; of all these the data pertaining 
to "high” A/V ratio of 11.7 sec ^ or more comes from motions recorded 
within 27 km of the epicentre. 

Table 3.4 also shows that A/V ratio is generally higher for the 
high PGA group data than for the medium or low PGA group data; i.e., the 
expected trend of decrease in A/V ratio with increase in the epicentral 
distance is generally seen. The only exception to this are the 
horizontal ground motion records of event "Oct" and the vertical ground 
motion records of event "April". However, the variation in A/V ratio 



with distance is not so significant for the Indian data as it is for the 
other earthquakes. For instance, in the present case, the epicentral 
distance ranges frcxn about 5 km to 320 km; for such distances, A/V ratio 
for earthquakes elsewhere could vary by more than 100%, which is not the 
case in the Himalayan data. 

Generally, the A/V ratio for vertical connpcaient of ground motion is 
comparable to that for the horizcmtal components. 

On the whole, a value of around 20 sec ^ could be considered as an 
appropriate mean value and 12 sec”^ as the 
mean-minus-one-standard-deviatioin value for the Himalayan earthquake 
records. Ihis is much higher than the values reccmmended in the 
literature based cn records elsetdiere. Ihis has very serious 
implications for earthquake-resistant design of important structures in 
the Himalayas. This means that the earthquake ground moticn in the 
Indian Himalayas are associated with (a) unusiaally low valxie of ground 
velocity and hence are far less damaging, (b) unxjsually high value of 
ground acceleration and hence are more damaging (particularly to the 
stiffer structures), or (c) a combination of the possibilities (a) and 
(b) above. To resolve this issue, systematic studies on attenuation 
relationships for the Himalayan earthqucJces are required for which far 
more authentic information on the earthquake source parameters is 
needed. In the meanwhile, for seismdc risk studies in India, a great 
deal of caution is required while using the empirical relationships 
(relating PGA and/or PGV with the earthquake magnitude and ^icentral 
distance) developed for data from elsewhere. 

Table 3.5 shows the mean value of A/V ratio obtained in each of the 
six Himalayan earthquakes. The results are tabulated in increasing order 
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of earthquake iragnitude. It is readily apparent from the table that mean 
A/V value decreases with the increase in magnitude. This is probably 
because, for a given PGA, the value PGV increases with the increase in 
earthquake magnitude. For the three earthquake magnitude ranges; (1) 
small magnitude (M < 5.5), (2) moderate magnitude (5.5 s M s 6.5), and 
(3) large magnitude (M > 6.5), the mean A/V value are also presented in 
the same table. 

3.3.2 AD/V^ Ratio 

Table 3.4 also shows the values of mean AD/V^ for the three groups 

in each event. 2U5/V^ for vertical component is higher than it is for the 

horizontal componaat; this is consistent with the data analysed by 

Mohraz (1976) and seen in Table 3.1. On the whole, the mean value of 
2 

AD/V ratio for the horizontal component is 4.6; this is lower than the 
value of around 7.0 for rock sites in Table 3.1, but higher than 2.3 to 
3.8 based on the Loroa Prieta data (Table 3.2). The data for "l^ril” has 
one record with unusually high value of AD/V^ (equal to 85.6) which 
could be due to a typographical error in the paper by Chandrasekaran and 

(1990). This gives unusually high values of the mean and 
mean-plus-one-standard-deviation for low PGA group of this event. Hence, 
in arriving at the overall mean and mean-plus-one-standard-deviation 
values for the vertical components (5.2 and 8.8, respectively), this 
single value has been disregarded. The data of Table 3.4 show no clear 
trend in the variation of AD/V with distance, except for the event 
"Aug". 


3.3.3 Ratio 
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A statistical sunnmary of the ratio of in the three PGA 

groups is presented in Table 3.4. This table shows that, except the 
event "Oct", the mean value of ratio is around 0.5 and 

mean-plus-ane-standard-deviation value is around 0.70. This event shews 
a rather strong vertical conponent of groxaid motion with a mean value of 
Ay/Ag as 0.66 and mean-plus-one-standard-deviation value of 0.9. For the 
entire data, the mean value of is 0.51 and the 

mean-plus-one-standard-deviation value is 0.70. These values are similar 
to those reported by Mohraz (1976) for data from elsewhere (Table 3.3). 
Thus, to be -on the conservative side, the design for vertical 

acceleration should be based on Ay/Ag ratio as around 0.67 or 0.75 as 
against the value of 0.5 presently being used in the Indian code. 

3.4 DORATION OF STRONG SHRKING 
3.4.1 Definition 

The total duration of a complete accelerogram represents the 
dxiration diuring which a earthquake ground acceleration exceeded a 
threshold or critical acceleration value set for the instrument. 

However, this duration may not represent the duration dviring which the 
structure was subjected to severe shaking. In the past, several 

procedures have been suggested for extracting a strong motion segment 
from an accelerogram; each of these give different value of the duration 
as seen in Fig. 3.3. As yet, no universally accepted definition for 
duration of strong motion is available. 

Bolt (1969) and Page et al . (1972) have suggested a "bracketed 
duration"; this is the duration between the first and the last 
acceleration peaks exceeding a specified value (usually 0.05g). Since 
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PGA decreases with increase in distance from epicenter, tfiis definition 
leads to a decrease in the duration of strong motion with an increase in 
distance from the source. 

Trifunac and Brady (1975) define the duration of strong shaking as 
the time duration during which the middle 90 percent (i.e., 5% to 95%) 
contribution of the "acceleration intensity" (defined as jA^dt) takes 
place. They observe that the average duration for soil is approximately 
10-12 sec longer than that for rock, and that it increases at the rate 
of 1.0 - 1.5 sec for every 10 km increase in distance from the 
shock-source. 

McCann and Shah (1979) define the duration of strong ground motion 
based on the concept that the strong motion part can be defined as the 
period during which the RMS, or average rate of energy, is highest 
relative to the rest of the record. In this method, first the root mean 
square (RMS) value is evaliwtod at each of tho dioitir.o<l j^oints in thr> 
accelerogram, this is referred to as cumulative root mean square 
function (CRF) . Rate of change of CRF during total duration of the 
record is then examined to give two cut-off points. The final cut-off 
time is obtained at a time where rate of change of CRF becomes negative 
and remains so for remainder of the record. The initial cut-off time is 
obtained in the same manner, except that the search is performed in the 
reverse direction of time scale starting from the end of the record and 
continued till beginning of the record. 

3.4.2 Comparison of Different Definitions 

A graphical comparison of duration as pei the Trifunac and Brady 
(1975) and McCann and Shall (1979) definitions is shown in Fig. 3.4 for 
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sore records from the Californian earthquakes. This figure shows that 
duration by the McCann and Shah method is shorter than that by the 
definition of Trifunac and Brady (1975). Cn the other hand, Table 3.7 
shows the average duration of records frcm the Lcma Prieta earthquake 
based cai three different definitions. This table indicates that the 
definition by McCann and Shah gives duration longer than that by the 
Trifunac and Brady method. This taJDle also shows that the duration is 
generally higher for alluvion sites than for the rock sites. 

3.4.3 Variation' in Duration with Distance and with R/V Eatio 

Tso et al . (1992) demonstrate that the records with high A/V ratio 
have shorter durations of strong shaking than those with intermediate 
A/V ratio, whereas the records with low fl/V ratio have longer durations. 
They explain this in terms of epi central distance as, "Larger the 
magnitude of an earthquake, the longer the duration of strong ground 
motion if the distance from the epicentre remains constant. If the 
magnitude of the earthquake is kept constant, duration increases with an 
increase in epicentral distance. Since the records with high A/V ratios 
were generated in the proximity of small or moderate earthquakes, 
whereas those with low A/V ratios were obtained at large distances frcm 
large or moderate earthquakes, the correlation of the A/V ratio with the 
duration of strong shaking can be expected.” 


3.5 DURATION OP STRONG SHAKING FOR IKE DATA 

The duration of strong shaking has been calculated for each of 
thel98 horizontal records based on the Trifunac and Brady (1975) 
definition (Appendix A). The relatiojaship between the duration and A/V 
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ratio for each of the events is shown in Fig. 3.5. The points for 
records belonging to different PGA groups (low, medium, and high) are 
shown with a different notation. The plots of Fig. 3.5 show no definite 
relationship or trond between duration and A/V or between duration and 
PGA group (i.e., epicentral distance), though such relationships have 
been reported in literature for studies based on data from other 
earthquakes (Tso et al . 1992) . 

This can also be seen in Table 3.8 which shows the mean duration of 
strong shaking in the three PGA groi?>s for the six Himalayan 
earthquakes. For "April" and "Fdb", the trend of increasing duration 
with decrease in PGA (i.e., increase in distance from shock-source) is 
quite clear. But, exactly opposite trend is observed in the event "Aug", 
where high PGA group records have greatest duratio of strong shaking. 
For events "Sept" and "Oct", the Table shows the largest duration of 
shaking in the medixm PGA group records. However, Fig. 3.5 clearly 
confirms the fact that larger the earthquake magnitiade, grater is 
duration of of strong shaking. For "sept" (M = 5.2) and "April" (M = 
5.5), duration of most of records is less than 10 sec, while for events 
"Oct" (M = 6.6) and "Aug" ( M = 7.2), duration of most of the records is 


more 10 sec. 
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Table 3.1 Sumnary of V/A and AD/V ratios for earthquakes in 
California (Mohraz 1976) . 


Site category 

Group 

Mean 

V/A 

(in/s)/g 

Approx, 
mean A/V 

(Sec”^) 

AD/V^ 1 

'Mean 

84.1 

percentile 

Rock 


27.0 

14.5 

6.9 

11.0 

S 

30.0 

13.0 

7.0 

11.2 


31.0 

12.6 

7.6 

11.8 

Less than 30 ft of 
alluvium underlain 
by rock 

n 

37.0 

10.6 

5.2 

7.7 

s 

43.0 

9.1 

5.2 

8.2 

V 

37.0 

10.6 

8.5 

13.3 

30 - 200 ft of 
alluviun underlain 
by rock 

L 

33.0 

11.8 

5.6 

7.8 

S 

41.0 

9.5 

4.3 

6.4 

V 

33.0 

11.8 

9.1 

13.7 


B 

51.0 

Wum 

4.3 

6.0 

Alluvivm 

s 

62.0 

6.3 

B 

4.4 


B 

51.0 


B 

BH 


L : Horizontal ccxiponents with the larger peak ground acceleration 
S : Horizontal carponents with the smaller peak ground acceleration 
V : Vertical coiponents 
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Table 3.2 Summary of mean V/A and AD/V^ ratios for different 
source to site distances in Loma Prieta earthquake 
of 1989 (Mohraz 1991). 



R : Distance of recording station from ruptured surface 


Table 3.3 Summary of the ratio of peak vertical to larger 
peak horizontal acceleration (Mohraz 1976). 


Site category 

Mean 

84.1 

Percentile 

Rock 

0.52 

0.69 

Less than 30 ft of alluvium underlain 
by rock 

0.49 

0.62 

30 - 200 ft of alluvixjni underlain 
by rock 

0.46 

0.66 

Alluvium 

0.45 

0.61 
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Table 3.4 Surmary of VV and AD/V^ ratios for the six <artbquake 
events in Himalayas. 


Event 

PGA 

Grotq? 

Horizontal 

ccxipments 


Vertical 

coqponent 


AD/V^ 

PQA 

Grot®) 

A/V (sec ^) 

adA^ 

m 

m - c 

m 

m + c 


m 

m - 

in 

m + <7 


HBBi 

21.1 

14.4 

2.9 

4.2 

H (5) 

23.2 

17.9 

5.6 

.7.8 

April 

M (4) 

21.4 

15.7 

4.2 

IQQI 

M (2) 

14.3 

13.8 

3.0 

3.5 


L (4) 

15.7 

7.2 

4.2 

6.7 

L (2) 

22.7 

7.0 

44.8 

89.5 


H (6) 

32.5 

25.8 

3.5 

4.4 

H (3) 

36.0 

30.1 

9.2 

13.6 

Sept 

M (8) 

31.0 

'■WMI 

5.9 

10.8 

M (4) 

26.8 

22.8 

10.1 

12.5 


LdO) 

26.5 

20.1 


8.6 

L (5) 

22.9 

1 

13.6 

6.9 

IB 


H (0) 


- - 


r- -- 

IQQI 



Bl 

— 

May 


23.7 

18.3 

6.5 

9.9 

M (6) 

19.5 

11.8 

?.3 

11.2 


1.(16) 

17.3 

12.6 

5.0 

6.5 

L (8) 

15.2 

12.5 

HQ 

6.5 


H (2) 

22.9 

21.3 

5.0 

■D 

H (1) 


— - 

9.2 

— 

Feb 

M(12) 

22.6 

17.4 

■g 

m 

M (6) 

25.1 

17.7 

8.1 

11.6 


L(22) 

18.7 

12.7 

5.2 

8.1 

L(ll) 

15.2 

9.2 

4.6 

6.2 


H(28) 

23.3 

17.1 

4.0 

5.9 

H(14) 

25.8 

18.5 

6.8 

9.4 

Aug 

M(32) 

16.7 

9.8 

3.6 

5.6 

M(16) 

16.3 

7.8 


7.5 


L (6) 

13.3 

9.4 

2.6 
t 

3.5 

L (3) 

12.8 

8.9 

3.2 

3.4 


H (8) 

14.3 

12.2 

Q 


H (6) 

15.0 

10.3 

2.9 

3.6 

Oct 

M(12) 

19.1 

13.0 

IB 

6.8 

M (4) 

14.4 

7.4 

2.9 

4.3 


h (6) 

14.7 

11.3 

3.2 

4.6 

L (3) 

10.7 

9.0- 

2.9 

i 3.2 

All 

(198) 

20.2 

12.5 

4.6 

8.0 

(99) 

19.6 

10.8 

5.2 

8.8 


* : Nunrber of records 

: High PGA grotp (PGA > 0.1 g) 

: Mediim PGA grov 5 > (0.1 > PGA > 0.05 g) 
; Low PGA grot®) (PGA < 0.05 g) 

: Mean 

+ 0 ; Mean-plus-one-standard-deviation - 
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Table 3.6 Sunmary of peak vertical to larger horizontal 
acceleration ratio for the data used in the 
present study. 


Event 

PGR 

Group 

V*H 

V*H = 

>f each event 


m 

m + 0 

m 

m + c 


H (5)* 

0.34 

0.46 

■1 


April 

M (2) 

0.57 

0.76 


0.74 


L (2) 

0.81 

0.97 

■1 



H (3) 

0.45 

0.47 



Sept 

M (4) 

0.41 

0.52 

0.44 

0.51 


L (5) 

0.40 

0.60 




H (0) 







May 

M (6) 

0.39 

0.52 

0.51 

0.72 


L (8) 


0.81 




m 

D3 

0.89 



Fdb 

m 

0.41 

0.47 

0.50 

0.72 


Uil) 

0.52 

0.77 



■ 

H(14) 

Rg 

! 0.55 

■■ 



M(16) 

0.49 

0.63 


0.61 

■ 

L (3) 

0.59 

0.67 

Bi 



m 

0.85 

1.05 



Oct 

M (6) 


0.77 

0.65 

0.90 


L (3) 

0.53 

0.80 



All 

(99) 



0.51 

0.70 


* : Number of records 

: Hig^ PCa groi^ (PGR > 0.1 g) 

: Medium PC3R group (0.1 > PGR > 0.05 g) 
: Low PGR group (PGR < 0.05 g) 

: Mean 

+ o : Mean-plus-one-standard-deviation 
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Table 3.7 Average duration of strong motion (in sec) for accelerograms 
recorded in the Lama Prieta earthquake of 1989 (Mohraz 1991). 


Category 

Page et al. 

Trifunac and Brady 

McCann and Shah 

Rock 

Alluvium 

Rock 

Alluvium 

Rock 

Alluvium 

Near-field 
(R < 20 km) 

11.6 

14.4 

8.5 

10.5 

10.0 

11.6 

Mid-field 

(20 km > R > 50 km) 

6.2 

12.5 

15.6 

18.6 

15.9 

19.8 

Par-field 
(R > 50 km) 

1.9 

6.7 

11.6 

11.5 

14.7 

14.8 
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B 

M 

L 


Table 3.8 Sumary of duration of strong 
shaking for the records obtained 
from the six Himalayan earthquakes. 



: Nunnber of records 
: High PGR group (PGR > 0.1 g) 

: Mediixn PGR group (0.1 > PGR > 0.05 g) 
: Low PGR group (PGR < 0.05 g) 
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EPICENTRAL DISTANCE (KMJ 

Fig. 3.1 Relaticjnship between earthquake magnitude (M) and 
epi central distance (R) for low, intermediate, and 
high A/V records (Tso et al. 1992). 



10 *’ 10 "' 10 * 10 ' 


PERIOD (SECOND! 


Mean 5% damped elastic acceleration response spectra 
for the three A/V groups and whole ensenfcle of 
records scaled to a coranan PGA of Ig 
(Tso et al. 1992). 


Fig. 3.2 
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60 
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Conpariscjn of strcang motion duration for the S69E 
component of the Taft, earthquake of July 21, 1982 
using different procedures 
(Mciiraz and Elghadamsi 1989). 


s 

-H 

4J 

«a 



Pig. 3.4 


Duration by 

Trifunac and Brady (1975) 

Ccnpariscn of duration of strong shaking from 
Trifunac-Brady (1975) and McCann-Shah (1979) 
definitions (McCann and Shah 1979). 



CHAPTER IV 

ELASTIC RESPONSE SPECTRUM 


4.1 nmiODOCTICK 

Response spectrum is a very useful way to characterize strung 
ground motion from the view point of engineering applications. Response 
spectnm of a ground motion is a plot showing maximum response induced 
that ground motion in single degree of freedom (SDC^) systems of 
different natural frequency (or nattiral period), but for a given value 
of damping. The response here can be any quantity of interest siich as 
the displacement, velocity, or acceleration. For engineering purposes, 
the following response quantities are most conmonly plotted (Fig. 4.1): 

(a) Maximum relative displacement (displacement spectrum) 

SD(«,C) = |u(t)|^ (4.1) 

t 

(b) Maximum relative velocity (velocity spectrum) 

SV(w,C) = l^(t)lmax 

t 

(c) Mautimum absolute acceleration (acceleration spectrum) 

SA(w,C) = |u(t) + (4.3) 

t 

(d) Pseudo spectral velocity (pseudo-velocity spectrum) 

PSV(«,C) = « SD(«,C) ' (4.4) 

(e) Pseudo spectral acceleration (pseudo-acceleration spectrum) 

PSA(M,C) = SD(w,C) 

where u(t), ii(t), ii(t) ace relative displacement, velocity, and 
acceleration, respectively, at instant t for a SDQF oscillator of 
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natural frequency « and dasiping C; and u^(t) is the groiaid acceleration. 

A response spectnm r^resents the coirbined influence of anplittde 
of ground motion, frequency content, and to some extent the duration of 
ground shaking, on different structiores (Seed and Idriss 1982). Most 
importantly, a response spectrum provides a most convenient way of 
evalxiating the maxinun acceleration, and thus neudnun inertia force, 
developed in a structure subjected to the given base motion, 
irrespective of vdiether structvure brfiaves as a SDOF system or MX®" 
system. Also, since a SDOF system responds mostly to the energy in the 
ground motion near its own natural frequency, the response spectxnxn 
characterizes the energy in ground motion. 

4.1.1 Tripartite R^resentation of Response Spectrum 

It is obvious from Egs. (4.4) and (4.5) that SD, PSV, and PSA are 
directly related and provide essentially the same information; hence the 
corresponding response spectra can be presented on a tripartite or a 
four-way logarithmic plot by a single curve for each value of damping. A 
typical response spectrun on a tripartite-plot is shown in Pig. 4.2, 
which is essentially a plot of PSV versus T on log-log scale with two 
additional log scales at 45 ** to period scale for reading SD and PSA. The 
parameters PSA and PSV have certain characteristics that are of 
practical interest. The PSV is close to SV for high frequencies, and 
practically equal for intermediate frequencies, but different for low 
frequencies. PSA is practically equal to SA for damping levels 
encountered in most ^ineering applications. 


4.1.2 Design Spectnzn 
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Response spectrxm, being a plot of naxinun response (of a SDCF 
oscillator with varying natural frequency) to a given ground motion 
record, is a description of that particular ground motion. Hence, its 
use is suitable for analysis of a structure subjected to that ground 
motion; but it may not be suitable for the purposes of desi^a (e.g. , 
Anderson 1989). Purthennore, a single earthquake record has a particular 
frequency content which gives rise to the jagged, sawtooth appearance of 
peaks and valleys as shown in Fig. 4.2. When such a spectrum is used for 
design, structuures with slightly different natural periods may reflect 
vastly different maximum response and hence the design parameters. To 
overcome this problem/ the concept of smooth-shaped response spectrum 
has been introduced for design. 

A usual method to obtain a smooth-shaped spectrum is through 
normalization and averaging of several response spectra obtained from 
records with conroon characteristics. Amongst the sugrgested normalization 
procedures for response spectrum, the most pqpular is normalization to 
unit peak ground acceleration; at times the normalisation is done with 
respect to other parameters such as peak ground velocity or 
displacement. Principles of statistics are then used to create a 
smooth-shaped spectrum. The mean or median shape of response spectra so 
obtained is generally used for earthquake-resistant design of normal 
building structures. Use of such a spectrum implies that for a given 
value of peak ground acceleration, there is about 50% probability that 
the maximum acceleration of the SDOF oscillator will be exceeded. High 
risk structures may be designed to higher level such as the 
mean-plus-one-standard-deviation shape, which implies that the 
probability of exceedance is only about 16%. To account for the effect 
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of nagnitude and distance of different likely events in the vicinity of 
a major project, site-specific design spectnm is sometimes dbtained by 

t 

drawing an envelop of more than one spectra (e.g.. Pig. 4.3). 

4.1.3 Hie Actual Design ^pectrun 

The above discussion pertains to the shape of design spectrun for a 
given normalised acceleration, velocity, or displacement. The acttial 
design spectrum is obtained by multiplying the entire spectrum shape 
with the design ground motion parameter, for instance the maxinutn ground 
acceleration (zero period acceleration, ZPA). The magnitude of this 
design parameter has to be decided considering the expected ground 
motion at the site, the factors of safety and load factors to be used 
in design, the ductility to be provided in the structure, the damping 
characteristics, and the procedure to be used for estimation of natural 
period for the stnicture (e.g., Housner and Jennings 1982). 

4.1.4 The C^jective of the Sttidy 

A statistical study of the response spectrum shape has been carried 
out for the records obtained from the Himalayan earthquakes, with a view 
to see if there are any appreciable differences in the ground motion 
characteristics of the Himalayan earthquakes vis-a-vis earthquakes 
elsevdxere. A comparison has also been made with the spectnm shapes 
currently in design use. 

4.2 REVIEW OP raST STUDIES OW DESIGK SPECTRDM 
4.2.1 Housner's Spectrum 

While Biot (1941) initiated the application of response spectnm 
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to earthquake engineering in 1941, it was not until 1959 that Housner 
developed the first smooth design spectra (Housner 1959). These (Fig. 
4.4) were obtained using eight strong motion records from four 
earthquakes recorded in California. The curves were arbitrarily scaled 
to 0.125g at zero period. Tlie Indian code (13:1893-1984) has adopted 
the shape given by Housner with minor modifications (Fig 4.5). 

4.2.2 Newroark and Hall S^pectrun 

Another set of design spectra, used extensively in the nuclear 
power plant design in late sixties and early seventies,, were suggested 
by Newmark and Hall (1969). They proposed a procedure for estimating 
site-independent design spectnxn, which was based on the fact that 
response spectrum over certain frequency ranges is related by an 
amplification factor to the peak value of ground acceleration, velocity 
and displacement. To appreciate this, consider the response spectrum of 
Fig. 4.2 and compare it with the peak ground displacement, peak ground 
velocity, and peak ground acceleration valties marked on the same plot. 
It is seen that at small frequencies, the maximum relative displacement 
is large, whereas pseudo-acceleration is very anall. it large 
frequencies, the relative displacement is extremely small, whereas the 
pseudo-acceleration is relatively large. At intennediate freguencies, 
the pseudo-velocity is stibstantially larger than at either eni of the 
spectrum. For this reason, three regions are ustially identified in a 
response spectrum as low-frequency or displacennent regim, the 
intermediate-frequency or velocity region, and the high-freqaency or 
acceleration region. In each region, the corresponding ground nation is 
attplified the most. Fig. 4.2 also shows that at small frequencies (0.05 
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Hz or less) the spectral displacement approaches the peak ground 
displacement, indicating that for very flexible systems the maximum 
relative displacement in a simple oscillator is equal to the ground 
displacement. At large frequencies (30 Hz or more), the 

pseudo-acceleration approaches the peak ground acceleration, indicating 
that for rigid systems the absolute acceleration of the mass is same as 
that of the ground. Newmrk and Hall suggested that the three 
amplifications (displacement, velocity, and acceleration) are constant 
in their respective regions and can be determined statistically. Tables 
4.1 and 4.2 give the amplification factors suggested Newnark and Hall 
for horizontal ground motions for median and 84 percentile cases. Thus, 
in this method, first the estimated values of peak ground acceleration, 
velocity, and displacement for the site are plotted on a tripartite 
logarithmic paper; then using the amplification factors, for the desired 
percentage of critical damping, the peak ground motion valiws are 
amplified to give a smooth design spe<^rum for the site. Their spectra 
for "standard” earthquake with peak ground acceleration, velocity and 
displacement of l.Og, 121.8 an/ sec and 91.4 csn, respectively, are shown 
in Fig. 4.6. It can be noted from the figure that curves for different 
damping values approach the peak ground acceleration at different 
period. Newmark and Hall observed that though the amplification factors 
for horizontal and vertical spectra remain unchanged, the vertical peak 
ground accelerations are about two-thirds of the horizontal peak ground 
accelerations. 


4.2.3 Study hy Blune and Mohraz 

After San Fernando earthquake of 1971, two independent studies, one 
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by Blume et al. (1972) and the other by Mohraiz et al. (1972), were 
carried out with the objective of developing recomnendations for 
horizontal and vertical design response spectra for nuclear facilities. 
These studies were based on a much larger number of earthquake 
grour»i-motion records than used ty the Newmark and Hall study. The 
results of these stvdies formed the basis for the guidelines contained 
in Nuclear Regulatory Coranissian (NRC) Regulatory Guide 1.6 (1973). The 
technique of r^resentaticai was very similar to the one proposed by 
Newmark and Hall, exc^t that these studies recognized that the 
anplifications in the three regions need not be constant always and may 
vary linearly. To define horizontal response spectra, the peak 
horizontal ground acceleration and displacement levels are first 
established. The peak horizontal ground displacement waus considered to 
be proportional to the peak horizontal grouni acceleration and was fixed 
at 91.4 am for an acceleration of Ig (for a standard earthquake). The 
bounds of each of spectrum are established by five line segments, 
similar to the Newmark and Hall procedure. The amplification factors and 
the response spectra for "standard" earthquake, for horizontal ground 
motion, as adopted by NRC are shown in Fig. 4.7. The vertical response 
spectra were constructed in a manner similar to horizontal response 
spectra, but three differences were incorporated into the procedure: (1) 
the frequency control points were shifted slightly, (2) the 

amplification factors were different, and (3) the value of peak 
horizontal acceleration is used as the initial reference valxre. The 
amplification factors and the respcmse spectra for vertical grotnd 
motion are shown in Fig. 4.8. The NRC spectra was also used extensively 
for seismic design of many non-nuclear facilities. 
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4.2.4 ^pectnm for Different Soil CScaiditions 

In 1976, two stiadies, one by Seed et al. (1976) and the other by 
Mohraz (1976), considered the influence of soil conditions on the shape 
of response spectnjm. The study by Seed et al. used a total of 104 
horizontal components of earthquake recoi:xis with peak ground 
acceleration greater than O.OSg from 23 different earthquakes. *flie 
records were divided into four categories based on soil condition at 
recording station. The mean and mean-plus-one-standard-deviatioo spectra 
for the four soil types, obtained by Seed et al. are seen in Figs. 4.9 
and 4.10, respectively. The figiires show that for periods greater than 
about 0.4 sec, the amplifications for rock are substantially lower than 
those for soft to mediim clay and for deep cohesionless soil. Ihe ABRB 
safety guide (S-11 1990) recommends standard response spectra for rock 
and soil sites based on the work of seed et al. (1976); Pig. 4.11 shows 
the shape recommended by AERB for rock sites. 

4.2.5 Study by Mohraz 

Mohraz (1976) considered a total of 162 records and divided these 
into four soil categories. To exclude the influence of the horizontal 
component of records with smaller peak ground acceleration on 
amplification, he divided the response spectra for each soil category 
into three sets: horizontal components with the larger of the two peak 
ground accelerations, horizontal components with the smaller of the two 
peak ground accelerations, and vertical components. Figure 4.12 shows 
the mean acceleration spectra for different soil conditions for the 
horizontal components with larger of the two peak ground accelerations. 
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is seen that the arnplificaticxi for alluvium extends over a larger 
freguency regicai than the amplifications for the other three soil 
categories. Also, for short periods, the spectral ordinates for alluviun 
lower than for the other soil types, whereas for intermediate and 
long periods, these are higher. 

Tlius, local soil conditions at the site significantly alter the 
freguency content in the ground motion and hence the response ^ectra. 
In the recent years, this has lead to the adqption of different design 
spectra for different soil conditions. Fig. 4.13 shows one such ^)ectra 
adopted by UBC (1991) . 

Mohraz (1978) studied the influence of earthguake nagnittn^e on 
acceleration amplification for alluvium. As seen in Fig. 4.14, 
acceleration amplifications are larger in earthquake of soagnitude 
jjetween 6 and 7 than for those with magnitude betweai 5 and 6. 

Mohraz (1991) examined the influence of source-to-distance on 
spectral shapes for the records obtained from the Loma Prieta earthjuake 
of 1989. The accelerograms obtained from rock and alluvium were divided 
into three subgroups based on distance of the recording station from the 
ground projection of the ruptured surface. The groups were labeled as 
’’near- field" for records with distance less than 20 km, "imid-field" for 
records with distance between 20 and 50 km, and "far-field" for records 
with distance greater than 50 km. Fig. 4.15(a) shows the spectral shapes 
for the three subgroups obtadned for 5% damping using horizontal 
coiTponents of accelerograms recorded on rock sites. The figure shows 
that the amplification for the near-field is substantially smaller than 
the anplification for mid- or far-field for periods greater than about 
0.5 sec. For periods less than 0.5 sec, however, the amplification for 
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the near-field is greater. Pig 4.15(b) shows the results of similar 
exercise for records obtained from alluvium. The figure indicates that, 
though distance influences the spectral shapes, the influence is not as 
prcaiounced as it is for rock. 

4.3 METHOD OF ANALYSIS 

The data set ccmsisted of 198 horizontal and 99 vertical cczrponent 
records, obtained from the six events referred to in chapter II. The 
ordinates of the pseudo-acceleration response spectra (PSA) are conputed 
by the Newmark-beta numerical method of solution (e.g. , Bathe and Wilson 
1978). The ordinates are evaluated for a period range of 0.0075 sec to 
3.0 sec at an interval of 0.0075 sec, thus generating about 400 points 
in this range. The conputations are carried out for 0, 0.5, 2, 5, 10, 
and 20 percent of critical danping. Since the peak ground acceleration 
for various accelerograms differ, the conputed responses cannot be used 
in statistical study without normalization. In this sttidy, normalization 
of spectra by peak ground acceleration is adopted. Thvjs, each of the 
pseudo-acceleration spectra is normalized by respective PGA, giving zero 
period ordinates of unity. 

The normalized spectra in each of the events are separated for 
horizontal and vertical components. These are then averaged to obtain 
mean and mean-plus-one-standard-deviation spectral shape for horizontal 
and vertical conponents representative for individual events. Similarly, 
a mean and mean-plus-one-standard-deviation spectra for horizontal and 
vertical component, considering all the records from 6 events is also 
evaluated. The average spectra are smoothened by adopting a 
five-point -averaging scheme. The "smooth-shapes" thus obtained are then 
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canpared with the spectral shapes curr^tly in lise. 

4.4 RESDLTS 

The mean pseudo-acceleraticai spectra for 6 events for dairping 
values of 0, 0.5, 2, 5, 10, and 20 percent of critical danping for 
horizontal and vertical conponents, are shown in Pig. 4.16. 

4.4.1 General shape 

A visual caiparison of spectra in the Fig, 4.16 suggests that in 
all the events,' for a given damping value, there is a build ip in 
acceleration amplification from unity at zero period to a peak value, 
which usually occurs in the period range of 0.1 to 0.3 sec. For periods 
longer than 0.3 sec the acceleration amplification drops down. However, 
for periods greater than about 2.0 sec the acceleration amplification is 
quite low and remains almost constant. 

For low damping values (less than 5%), the spectral shape and the 
value of peak amplification in different events vary significantly. 
However, for higher damping values, the peaks of spectra in different 
events are nearly same and the variation in the shape is also 
comparatively less. For instance, in the event "Sept" (M = 5.2) the peak 
amplification for 0% and 5% damping are about 4.8 and 2.4, respectively, 
while in "Aug" (M = 7.2) the corresponding amplifications are about 9.4 
and 2.6, respectively. Thus, a significant variation to the extent of 
200%, is observed in zero damping spectra while the difference is 
insignificant for 5% damping. Though the shape of spectra differs from 
event to event, the peak acceleration amplification in the six events 
for 5, 10, and 20 percent of critical damping is nearly same and is 
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about 2.6, 2.0, and 1.4, respectively. 

4.4.2 Variaticn within the Data 

l^pectra in the "Sept" event drops down suddenly at period of about 
0.3 sec. This could be due to its low magnitude, ’^y" and "Feb" are of 
almost same magnitude, their spectral shapes are in general similar, 
except for some variation in maximum amplification for danping less than 
5% of critical. "Aug" is of the largest magnitiKle, and in this 
earthquake most stations were at large epi central distances (greater 
than 150 km). The spectra for this event show corparatively (with 
respect to spectra for the other earthquakes) more energy in 
intermediate period range (0.4-1. 2 sec); this is expected because as 
distance from the shock-source increases, the high period waves 
dominate. In case of "April", though the epicentral distances involved 
were quite low (less than 30 km), the ordinates in intermediate period 
range show a substantial anplifi cation. This is unexpected, since 
generally seismic waves in near-field (low epicentral distance) have 
high frequency content. "Oct" spectra are unique in the sense that cxily 
these show a flat top as compared to sharp tips observed in spectral 
shapes of the other events, 

4.4.3 Variation in Shape of Spectrum with Earthquake Magnitude 

Considering the large variation in arrplifications for low values of 
danping, and the fact that in the design of ordinary buildings a danping 
value of 5% is usually adopted, the variation in spectral shape with 
earthquake magnitude has been studied for 5% danping. Fig. 4.17 shows 
the mean spectra for the six events for 5% danping. It is seen that 
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acceleration anplifications in the period range of 0.3 to 1.5 sec are 
highly dependent on earthquake magnitude. Larger the nagnitude, greater 
is the anplificaticn in this range. For instance, at natural period of 
0.8 sec, the acceleration anplificaticn is 0.3 for the event "Sept" (M = 
5.2) and it is 1.0 for the event "Aug" (M = 1 . 2 ), giving a variation of 
as much as 500%. For periods less than 0.2 sec the variation is 
negligible, while it is not very significant for periods longer than 1.5 
sec. This figure also suggests that one could obtain somewhat 
representative shapes of the response spectrvro based on earthquake 
magnitude. Fig. 4.18 shows shapes suggested by this data for rock-sites 
for (a) small -magnitude earthquakes (M < 5.5), (b) moderate-size 

earthquakes ( 5.5 s M i 6.5), and (cl large earthquakes (M > 6.5). These 
have been obtained from (a) the mean shape of event "Sept" (M = 5.2), 
(b) average of the mean shapes of evaits "April" (M = 5.5), (M = 

5.7), and "Feb" (M = 5.8), and (c) average of the mean shapes of events 
"Oct" (M = 6.6) and "Aug" (M = 7.2). 

The variation of acceleration anplificaticn in the intermediate 
period range with earthquake magnitude can also be anticipated from the 
A/V ratio. As seen in table 3.5, in the low-magnitude earthquakes, the 
value of A/V is high, i.e. , for the given value of PGA, the value of PGV 
is low. In the intermediate period range (the velocity region), the PGV 
governs the spectral ordinat^, and hence, low acceleration 
amplification are obtained in this range. 

Not many studies seem to have been done on the shape of 
pseudo-acceleration spectnan with earthquake magnitude. The study by 
Mohraz (1978) (Fig. 4.14) shows this variation but for alluvium sites. 
While the study by Mohraz based on data from elsewhere shows that ground 
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motions of large magnitude earthquakes have significantly larger values 
of acceleration amplification in all the three regions (acceleration^ 
velocity, and displacement), the data of the Himalayan earthquake (Fig. 
4.17) suggests that this is so only for the velocity region. 

It is obvious that when choosing spectnm shape for preparing 
site-specific design spectrum, the earthquake magnitude shoixld be 
explicitly considered, particularly in the intermediate period range; 
this is the period range in which a very large niiriber of industrial 
structinres lie. Presently, such a consideration seams to be absent in 
developing the design criteria for important projects in India. 

4.4.4 Hhat does Mean S^>ectra from Data R^resent ? 

Mean horizontal pseudo-acceleration spectra, derived by averaging 
all 198 horizontal componoit spectra is shown in Fig. 4.20. Comparison 
of shape of mean spectra from all horizcntal records with spectral 
shapes from "Oct" (M = 6.6) and "Feb" (M = 5.8), shown in Fig. 4.21, 
indicates that mean shape (from all horizontal records) lies in between 
the mean spectra for two events, suggesting that mean spectral shape 
from all records is a represoitative of earthquake of magnitude of about 
6.2 for moderate distances (80-130 km). Incidentally, the average of 
magnitudes of events studied is 6.0. 

4.5 VARIATION IN SPECTRAL SHAPES WITH DISTANCE 

To study the variation in response spectrum shape with epicentral 
distance, mean response spectra for high, medium and low horizontal PGA 
groups are plotted for different events in Fig. 4.19. The shapes for 
events "Sept", "May", and "Oct" do not show significant variation with 
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distance even though the peaks in the acceleration region vary somewhat. 
For "J^ril”, the ordinates at large values of natural period are higher 
for the low PGA grovp, but in this event epicentral distance for all the 
stations lie within 30 km. For event "F^”, the shape for high PGA group 
is significantly different from that for the other two grotps, but in 
this case only two records (i.e., one station) falls in the high PGA 
group and therefore this comparison may not be very general. In the 
intermediate period range, the events "Feb" and "Rog” show significant 
variation for different PGA grotips; the high PGA grov®) (low epicoatral 
distance) has much lower values of acceleraticHi amplification in this 
period range cis is also expected. 

Thus, in general, the Himalayan data does not show a significant 
variation in the spectral shape with epicentral distance; this is in 
contrast with the data from elsevdiere. For instance, compare the 
variation for rock sites from the Lena Prieta earthquake (Fig. 4.15) 
with the event "Oct" (Fig. 4.19). While the Loma Prieta earthquake shows 
a very significant variation even within the epicaitral distance of 80 
km, the data of "Oct" pertaining to the epicentral distances of 25-150 
km shows not much variation. Only for event "Aug", where large 
epicentral distances (100-320 km) are involved, there is significant 
difference in spectral shape with distance, that too only in the 
velocity region. 

On the whole, it could be ccaicluded that the Himalayan data shows 
somewhat insignificant variation in spectral shape with epicentral 
distance, except in the intermediate period range for the very large 
distances (say more than about 150 tan). This is in line with the 
observation made in Chapter III that A/V ratio for the Himalayan 
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earthquake does not vary as significantly with distance as the data from 
elsewhere would suggest. Therefore, it is obvious that the geological 
featxires in the Himalayas are such that the attenuation is not 
significantly different for waves of different frequencies. 

An examination of mean spectral shapes in the three PGA groL?>s 
together with mean A/V of the corresponding group (T^le 3.4) reveals 
that spectra with high A/V value have least amplification in moderate 
period range, whereas for spectra in low A/V groi 5 > just the opposite is 
true. This again demonstrates that A/V is a good indicator of frequency 
content of ground motion. 

4.6 MEAN-PLUS-CHE-STANDARD-DEVIATION SPBCTROM 

While the mean spectrum shape is often recommended for design of 
buildings, inportant facilities such as nuclear power plants aure usually 
designed on the basis of roean-plus-one-standard-deviation shape. This 
ensures that the variation in spectra ordinates for a given value of PGA 
is accounted for in an adequately conservative manner. In the case of 
Himalayan data, where many parameters such as local soil conditions, 
epicentral distance, and earthquake magnitude are not very reliably 
known, the mean-plus-one-standard-deviation shape has the added 
significance such a shape based on all the data accounts for variation 
in these factors. 

Fig. 4.22 shows the mean-plus-one-standard-deviaticn spectra for 
the six events. In general, shape of mean-plus-one-standard-deviation 
spectra is very similar to shape of mean spectra (Pig. 4.16). Siidlar to 
mean spectra for the three range of earthquake magnitude (Fig. 4.21), 
mean-plus-one-standard-deviation spectra are also obtained. Pig. 4.23 
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shows the average representative shapes of the spectra for large 
magnitvide (M > 6.5), moderate magnitude (5.5 ^ M ^ 6.5) and low 
magni tilde (M < 5.5) earthquakes. Fig. 4.23 reveals that the 

amplification in intermediate period range is highly depaident on 
earthquake magnitixie. As magnitude increases, the spectral ordinates 
also increases and the increase could be as high as 500%. 

Figure 4.24 shows the mean-plus-one-standard-deviation spectra 
evaluated considering all the horizontal ccnponoit records; these are 
the shapes that may be adopted for major projects in the Himalayas. 

Figure 4.25 shows the coefficient of variation (COV) of response 
spectra for 0, 5, and 20 percent damping, considering all the horizontal 
conponent records. The figure suggests that OOV decreases with increase 
in damping. Also, COV is highly dependent on the natural period. The 
value of OOV for 5% damping is as low as 0.17 at period of 0.1 sec and 
increases to a value of around 0.8 at the period 1.2 sec. This is 
because the scaling of spectra has been performed with respect to the 
PGA. 


4.7 COMPARISON OF SPECTRAL SHAPES 

The design spectrum adopted by IS:1893-1984 is compared with the 
the mean horizontal spectral shape for 5% damping in the six Himalayan 
earthquakes (Fig. 4.26). It is readily apparent from the figure that 
irrespective of magnitude of earthquake, the IS code spectra 
underestimates the acceleration amplification in period range of 0.1 to 
0.35 sec while overestimates in the period range greater than 0.35 sec. 
The mean spectra of event "Aug" (M = 7.2) is closest to IS spectra for 
periods greater than 0.5 sec. In this event, records were obtained from 
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large epicentral distances (100-320 km); this suggests that for long 
period range, the IS spectra actually represents large-magnitude, 
far-field spectra. IS code spectrum shape in the low-period range needs 
to be revised upwards in light of this data. will also be seen in the 
next chapter, this is the period range in which ductility is scmewhat 
less effective in reducing the seismic force on a structure. 

4.7.1 Mean Spectrum with the IS code and DBG (1991) ^pectnm 

A conparison of the mean spectral shape based on all 198 records 
with the design spectra of IS: 1893-1984 and UBC-1991 (rock site) hais 
been shown in Fig. 4,27. The shape given in UBC-1991 appears quite 
conservative as compared to the Himalayan data. Note however that 
UBC-1991 does not specify a value of damping for this spectruB shape; it 
is assumed here to be 5%. 

4.7.2 Mean-Plus-One-Standard-Deviaticai Spectra with the Spectra 

Fig. 4.28 shows cctrparison of mean-plus-one-standard-deviation 
spectral shapes from the data with the "standard spectra" suggested by 
AERB for rock sites. It is seen that except for zero dainping, the AERB 
spectra match very well with the mean-plus-one-standard-deviatian 
spectra obtained from the data. 

4.8 STUD? ON AMPLIFICATION FACTORS 

Representation of spectra in the characteristics trapeaaidal form 
as done in the previous works by Newmark and Hall (1969), Mokaz et al . 
(1972), Blume et al.(1972), and Mohraz (1976) has also been attempted in 
the present study. For this the ordinates of PSV , PSA, and spectra 
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are evaluated using Newmark-beta method (Bathe and Wilson 1978). The 
ordinates were ccnputed for a set of 200 periods, distributed between 
0.003 sec to about 30.0 sec. The intervals of periods are chosen such 
that they get tmiformly distributed on a log scale. 

From the log- log plots of mean and mean-plus-one-standard-deviation 
normalised PSA, PSV, and SD spectra in each of the events, the most 
sorplified region is identified. The period limits between vMch these 
amplifications occur are recognised as amplification regirai of that 
spectra. For instance, the most amplified region in a PSA spectra is 
recognised cus acceleration region. The spectrum in the most asrplified 
regicn is then represeited by a straight line, which in most cases is a 
good approximation. This representation is illustrated in Figs. 4.29 and 
4.30, which show mean PSA and PSV spectra, respectively, for "Oct" 
event. Also shown in those figures are the straight line representation 
of the spectra in the acceleration and velocity regions, respectively. 

Figure 4.30 also suggests that spectra in the displacement region 
(which is mostly between 3 to 10 sec) cannot be represented as a 
straight line as has been done in the previoias studies. In fact the 
shape of spectra in this region is unexpected. Further investigations 
are necessary for the representation of spectra in this regirai. Hence, 
in the present study anplification factors only for the acceleration and 
velocity regions are evaluated; this is sufficient for most of the 
engineering structures. 

The spectra in acceleration region are obtained as straight lines 
parallel to period axis, that is, constant anplification (constant with 
fn r.Arind') of peak oround acceleration is obtained. However, in 
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inclined straight lines (e.g. , Fig. 4.30), indicating that in this 
period range the airplification of PGV is linearly varying (on log-log 
plots) and it decreases as period increases. For each of the events 
considered, the mean acceleration and velocity anplificaticm factors are 
presented in Table 4.3 and Table 4.4. For velocity region, the spectra 
being linearly varying, amplification factors at the ends of the 
velocity regicm are presented. The I^les also indicate the period 
limits of the acceleration and velocity region as seen for these 
events . 

Table 4.3 -reveals that for low damping values (less than 5%) the 
acceleration amplification factor varies for the six events, however 
for higher damping values this factor is quite consistent for the events 
considered. Also, the period limits defining the acceleration region do 
ot vary significantly for different evaits. Generally, acceleration region 
of the Himalayan data is between 0.12 sec to 0.26 sec. Table 4.4 
indicates large variation in amplification factors in the period range 
defining the velocity region. Though the lower period end of velocity 
region is somewhat consistent in the six events, the higher period end 
varies significantly. This also gives rise to large variaticn in 
amplification factor at the higher period end of velocity region. The 
velocity region for the Himalayan data lies between 0.26 to 2.6 sec. 

Acceleration and velocity amplification factors evaluated 
considering all the horizontal component records, are compared with the 
amplification factors reported in previous stx:uiies by Newmark and Hall 
(1969), Mohraz (1976) and those recommended in NRC guide 1.6 (1973). The 
comparisons for 5% damping are shown in Pigs. 4.31 and 4.32. The 
comparison is also summarized in Tables 4.5 to 4.8 for 2, 5 ,10, and 20 



percent danping. 


Figure 4.31 for mean acceleration and velocity anplification 
factors, clearly reveals that mean acceleration airplification obtained 
from the Indian data is higher than that for data from earthquakes 
elsewhere. On the other hand, mean amplification in the velocity region 
is lower than those obtained by Newroark-Hall and Mohraz, Also the 
velocity amplification factor reduces sharply with period in the 
Himalayan ground motions. However, mean-plus-one-standard-deviation 
acceleration amplification from Indian data is comparable (Fig. 4.32) 
with those obtained by Newmark-Hal 1 , Mohraz and the one recommended by 
NRC guide. This means that the standard deviation in acceleration 
amplification in the acceleration region is much lower for the Himalayan 
data than for data from elsewhere. 

4, 9 SPBCTROM INTENSITT 

Peak ground acceleration is the most commonly used ground motion 
parameter in design. However, it has come under much criticism for its 
lack of correlaticxi with observed structural performance. Hence, many 
other parameters have been proposed to measure the intensity of 
earthquake ground shaking. One such parameter is spectrum intensity, 
proposed by Housner (1959), Spectrum intensity (SI) is defined as the 
area under a pseudo-velocity spectrum curve between the periods 0.1 to 
2,5 sec, i,e. 


„2.5 

SI O) = S, (T,^) dT 

'’O.l ^ 


(4.6) 


Where is PSV and P the percentage of critical danping. 



Tso et al. (1992) studied the correlatican between h/V ratio and 
spectral intensity, with ^ being 5%. They observed that mean spectral 
intensity evaluated from PSV curves normalized to a conmoe PGR of l.Og 
for low fl/V (A/V less than 7.9 sec group of records is over four 
times higher than that for the high A/V group (A/V greater than 11.8 
sec ^). Further, they show that SI correlates well with peak groxmd 
velocity (V), irrespective of A/V ratio (Fig. 4.33). They have obtained 
an aipirical relationship between SI and V as 

SI = 3.3V (4.7) 

For the data set used in this sttidy, a similar correlation between 
SI and V is obtained. Spectrum intensity for 5% damping, for all the 
horizontal component records used in the study are presented in 2^»pendix 
A for ready reference. Fig. 4.34 shows the plots of SI versus peak 
ground velocity, for the six earthquakes; also shown in the figure are 
the best-fit linear (in the form y = m x) relationships. Bnnpirical 
relationships for each of the six events are presented below 
SI = 2.65V (For "April", M = 5.5) 

SI = 2.19V (For "Sept”, M =5.2) 

SI = 3.36V (For "May", M =5.7) 

SI = 3.00V (For "Feb" M = 5.8) (4.8) 

SI = 3.56V (For "Aug", M = 7.2) 

SI = 3.98V (For "Oct", M =6.6) 

SI = 3.55V (For All records) 

The relationships indicate that, in general, slope of straight line 
correlating SI with V increases with increase in magnitude of 
earthquake. However, no definite relationship in slope and earthquake 
magnitude can be concluded from this study . 
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Table 4.1 Spectrvum amplification factors for horizontal 
elastic response (Newmark and Hall 1969). 


Damping , 

% Critical 

One Sigma (84.1%) 

Median (50%) 

A 

rv— 

D 

A 

V 


0.5 

5.10 


3.04 

3.68 

2.59 

2.01 

1 

4.38 

3.38 

2.73 

3.21 

2.31 


2 

3.66 

2.92 

2.42 

2.74 

2.03 

1.63 

3 

3.24 

2.64 

2.24 

msm 

'■KOI 

1.52 

5 

2.71 

2.30 

2.01 

2.12 

1.65 

1.39 

7 

2.36 


1.85 

1.89 

1.51 

1.29 

10 

1.99 

mmm 

1.69 

1.64 

1.37 

1.20 

20 

1.26 

IJLZJ 

1.38 

1.17 

msEm 

1.01 


A : Acceleration 
V : Velocity 
D : Displacement 


Table 4.2 Equations for Spectrum Amplification Factors 
for Horizontal Motion (Newmark amd Hall 1969). 


Quantity 

Cumulative 
Probability, % 

Equation 

Acceleration 

Velocity 

Displacement 

84.1 

- - - f - — 

4.38- 1.04 InO) 

3.38- 0.67 InO) 
2.73-0.45 InO) 

Acceleration 

Velocity 

Displacement 

50.0 

3.21-0.68 InO) 
2.31-0.41 InO) 
1.82-0.27 InO) 


p : Fraction of critical damping 
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Table 4.3 Sutmary of mean acceleratim anplification factors 
for the six Himalayan earthquakes. 



period 

limits 

Percentage of critical danping 


(sec) 

0.5 

2.0 

5.0 

10.0 

J^ril 

0.13-0.24 

4.2 

3.3 

2.6 

2.1 

S^t 

0.12-0.21 

4.1 

3.3 

2.6 

2.1 

May 

0.12-0.27 

5.8 

3.9 

2.8 

2.0 

Feb 

0.11-0.24 

4.7 

3.7 

2.6 

2.0 

Aug 

0.11-0.27 

5.4 

3.6 

2.5 

1.9 

Oct 

0.14-0.29 

5.4 

3.8 

2:1 

2.1 

All 

records 

0.12-0.26 

5.0 

3.5 

2.6 

2.0 


Table 4.4 Surmary of mean velocity anplification factors for 
the six Himalayan earthquakes. 



period 

limits 

Percentage of critical damping 

Event 

(sec) 

0.5 

2.0 

5.0 

10.0 

April 

0.24-0.90 

2. 8-2. 8 

2. 3-2. 3 

1.8-1. 8 

1.4-1. 4 

Sept 

0.21-0.30 

3. 6-3. 6 

2. 7-2. 7 

2. 3-2. 3 

1.6-1. 6 

May 

0.27-2.10 

3. 6-2.0 

2.7-1. 6 

1.9-1. 4 

1. 5-1.0 

Pd3 

0.24-2.00 

2. 9-1. 9 

2. 2-1. 6 

1.6-1. 4 

1. 3-1.0 

Aug 

0.27-1.20 

3. 9-3. 6 

1 

2. 7-2. 7 

2. 0-2.0 

1. 5-1.5 

Oct 

0.29-1.28 

3. 8-2. 6 

2. 9-2.1 

2.1-1. 7 

1. 6-1.3 

All 

records 

0.26-2.6 

3. 6-1. 7 

2. 6-1.3 

2. 0-1.1 

1. 4-0.9 




Table 4.5 Carparitive sinmary of mean acceleration anplification 
factors. 


Percent of 

critical 

damping 


Present INewmark & HalllMoBfaz (1976) 

for rode sites 


Period, Period, Period, 

0.12-0.26 sec 0.12-0.50 sec 0.12-0.33'Sec 


Table 4.6 Comparitive stitmary of mean velocity amplification 
factors . 


Present i Newmark & Hall) jHohraz (1976) 

for rock sites 


Percent of 

critical 

damping 


Period, 

0.26-2.60 sec I 0.50-3.33 sec 


2. 6-1. 3 


Period, 
0.33-3.33 sec 


2 . 0 - 1.1 


1.4-0. 9 


1 . 1 - 0. 9 
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Table 4.7 Camparitive sutmary of roean-plvis-one-standard-deviation 
acceleration 2ffrplification factors. 


Percent of 

critical 

damping 

Present 

study 

Newmark & Hall 
(1969) 

Mbhraz (1976) 
for rock sites 

NRC 

Period, 
0.12-0.28 sec 

Period, 
0.12-0.50 sec 

Period, 
0.12-0.33 sec 

Period, 
0.11-0.40 sec 

2 

4.6 

3.66 

3.80 

3.54-4.25 

5 

2.9 

2.71 

2.82 

2.61-3.13 

10 

! 2.4 

1.99 



2.11 

2.27-2.72 

20 

1.7 

1.26 

1.54 

1.90-2.28 


Table 4.8 Ccrnparitive sutmary of mean-pltxs-one-standard-deviation 
velocity arrplification factors. 


Percent of 

critical 

damping 

Present 

study 

Newmark & Hall 
(1969) 

Mohraz (1976) 
for rock sites 

Period, 
0.28-2.54 sec 

Period, 
0.50-3.33 sec 


2 

3. 6-1. 9 

2.92 

2.44 

5 

2. 6-1. 6 

2.30 

1.90 

10 

1. 9-1.1 

1.84 

1.48 

20 


1.37 

1.11 








































%>0ctral Velocity, In/sec 
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Fig. 4.1 SDOF System, siibjected to earthquake ground motion. 



Fig. 4.2 A typical response spectrum with a tripartite plot. 

Also, marked are the values of peak ground 
acceleration, peak ground velocity, and peak ground 
displacement (Gt^pta 1950). 


Velocity, Sk (M per sec) 
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I 
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Fig. 4.3 Idealized undarrped velocity spectnin ctirves 

illiastrating the effect of magnitude and distance, 
curve A, 25 miles from centre of large earthquake; 
ciurve B, 70 miles from centre of large shock; curve 
C, 8 miles from centre of small (M=5.3) shock 
(Housner 1970). 
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Fig. 



VELOCITY SPECTRUM CURVES 



anooth design spectr\m curves developed by Housner 
( 1959 ). 
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SPECTRAL VELOCirf, IN CENTIMETERS PER SECOND 
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PLRIOO.IN^tCONOS 



Fig. 4.6 "Standard" site independent horizontal response 

spectra modified fran Newmark and Hall 1969 (Ilayn 
1980). 




SPECTRAL VELOCITY. IN CENTIMETERS PER SECOND 


74 


PERIOD, m SECONDS 



Horizontal design spectra and relative values of 
spectrum anplification factors for control points, 
adopted by NRC regulatory guide 1.6 (1973) (Hays 


Fig, 4.7 




















SPECTRAL VELOCfTY. IN CENTIMETERS PER SECOND 
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PEfilOO. IN SECONDS 


20 1.0 0.5 0.2 0.1 0.05 0.02 



% of 

critical 

danping 

Acceleration 

Displacement 
D(0.25 Hz) 

A’ (50 Hz) 

A(33 Hz) 

B(9 Hz) 

C(3.5 Hz) 

0.5 

0.67 

1.0 

4.96 

5.67 

2.17 

2.0 

0.67 

1.0 

3.54 

4.05 

1.67 

5.0 

0.67 

1.0 

2.61 

2.98 

1.37 

7.0 

0.67 

1.0 

2.27 

2.59 

1.25 

10.0 

0.67 

1.6 

1.90 

2.17 

1.13 


Vertical design spectra and relative values of 
spectnm aitplification factors for control 
points, adopted by NRC regulatory guide 1.6 (1973) 
(Hays 1980). 


Fig. 4.8 
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Fig- 4.9 
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Fig. 4.11 flERB "standard" spectra for rock sites. 
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PERIOD, SEC. 

Average acceleration anplification for 2 percent 

damping (Mohraz 1976) . 
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Fig. 4.15 
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Fig. 4.17 Mean horizontal acceleration spectra for 5% danping 
for the six Himalayan earthquakes. 
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Mean horizontal acceleration spectra for 5% dating 
to tte three KR groups, for the six Hirnalayan 

Km lakes * 




Fig. 4.19 
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Fig. 4.21 
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Fig. 4.24 
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Mean-plus-one-standard-deviation horizontal spectra 

considering all records together for 0, 0.5, 2, 5, 
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Fig. 4.25 OOV of acceleration spectra for 0, 5, and 20% dairping 
considering all horizontal records together. 



Fig. 4.26 Corrparison of 13:1893-1984 spectra with mean spectra 
from six Himalayan earthquakes. 
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Coiparison of mean-plus-one-standard-deyiation 
anplification factors (a) acceleration anplification 
and (b) velocity amplification. 


Fig. 4.32 
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Fig. 4.33 


Correlatim between spectnm intensity and PGV for 
three groups of records having low, intermediate, and 
high A/V ratios. 














CHAPTERV 

DUCTILITY REDUCTION FACTO 


5.1 INIRODOCTICH 


Elastic response spectrum is obtained by assvndng that the 
structure behaves linearly elastic even during severe earthquake ground 
motion. However, it will be inpractical, or at least uneconomical, to 
design a structure to remain elastic during severe earthquake ground 
motion of rare occnrence. Hence, the seismic design specifications for 


structures must incorporate the inelastic response appropriately. It is 
well understood and recognized that when inelastic deformations are 
permitted in design by providing adequate ductility, the structure can 
be designed for a much lower seismic force than otherwise. 


The cod^, particularly in the earlier years, provided for seismic 
design spectrun which was much lower than the elastic spectrum for the 
expected severe ground action of rare occurence. Such an approach 
inplicitly assuned the inelastic response of structure. However, the 
current trend of seismic code develoEnent is to prescribe a realistic 
level of elastic design spectrun corresponding to severe ground motion. 
Force obtained by this spectrua is then reduced by a "response reduction 
factor" (R). 18:1893-1984 still has 'the first approach; many recent 
codes (Rrc-1978, OBC-1991, HEHRP-1991, etc.) have adopted the later 
approach. The response reducticc factor primarily consists of two 
factors, one factor is to reduce the elastic demand force to the level 


of maxinam 
depends on 
dissioation 


yield strength of the structure; this reduction basically 
ductility of the structure and hence on its energy 
capacity. Ibis factor is called the ductility reduction 
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factor, R^. The other factor, called overstrength, is due to the 
overstrength inheraitly introduced in the code-designed structures and 
is defined as the ratio betweai the maxirautn lateral strength of the 
structure and the code prescribed unfactored design base shear force. 

Fig. 5.1 shows idealized linearly elastic-perfectly plastic (EPP) 
global response of a building in terms of base-shear coeffecient versus 
roof displacement relaticsiship during severe ground shaking. Also shown 
in the figure is the response of building, if it retrains elastic. The 
naxinun base shear that develops in the structure, if it were to renain 
elastic, is where W is the total weight of building. During 
nonlinear response, the naximun roof displacemait of the structxire is 
'^nax* global ductility denand \i, defined as 



can be achieved through appropriate design and detailing of structure, 

then such a building needs to be designed only for maximum yield 

C 

strength of CyW, where Cy = ^ . Hence, plays an inportant role 

in evalxiating the design force on a structure. 

Ductility reduction factor (I»F, R^) depends on natural period and 
damping of the structure, its ductility, its load-displacement 
behaviour, and on the ground motion characteristics. In the present 
study, the structure has been assuned as a SDC®* elastic-perfectly 
plastic (EPP) system. The value of DRF has been taken as the ratio of 
maximum elastic response (PSA) to the maximum inelastic response (PSA) 
with ductility fi for a SDOF EPP oscillator with 5% damping, i.e. , 
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PSA (T, C = 5%)^^^ 
PSA (T, C = 5%, M) 


(5.2) 


In this chapter, the value of DRF has been studied for the 
Himalayan ground motion records, and compared with studies available in 
the literature for ground motions recorded elsewhere. 

5.2 REVIEW CF THE PAST STODIES 

The first attannpt to relate with fi was made by Newtark and Hall 
(1973, 1982) for SDOF systems with EPP force-deformation model. Based on 
analytical studies, they concluded that: 

- In structures of natural period less than 0.1 sec, ductility does not 
help in reducing the response. Hence, for such structures, no 
ductility reduction factor should be used. 

- For moderate period structures, corresponding to acceleration region 
of elastic response spectrum (T = 0.1 to 0.5 sec), the energy absorbed 
by an inelastic structure at its rnaxiiam displacement approximates 
that absorbed by an elastic system. This results in the ductility 
reduction factor as (Fig. 5. 2a) 

= (2/i - 1)^/^ (5.3) 

- For relatively long period structures, corresponding to velocity 
region (T = 0.5 to 3 sec) and displacemait region ( T = 3 to 10 sec) 
of elastic response spectrum, the maximum displacement of elastic and 
inelastic systems is nearly egml. This gives the value of ductility 
reduction factor as (Fig. 5. 2b) 


(5.4) 
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Figure 5.3 shows the DRP spectra proposed by Newrrark and Hall (1973). 

Riddell and Newrnark (1979, 1980) conducted a statistical analysis 
of the response of SDOF systems with elasto-plastic, bilinear, and 
stiffness degrading load-deformation model. For EPP, SDOF systems they 
obtained empirical fomulae for R^ depending on ductility (m) and 
viscous damping ratio (C) for different period ranges. They observed 
that: 

- The ordinates of the mean inelastic spectra do not vary significantly 
whoi various nonlinear models are used; differences occur nainly for 
intermediate frequencies and large ductility, and are practically 
negligible at the low and high frequsicy ends of the spectnin. 

- Use of the eleisto-plastic idealization provides, in almost every case, 
a conservative estimate of the average response to a number of 
earthquake motions. 

- It is particularly significant that, on the average, the stiffness 
degradation phenomenon is not as critical as one might expect a 
priori . 

Using ten "standard" ground motions from five different seianic 
events, Mahin and Bertero (1981) have computed displacement ductility 
demands for EPP SDOF systems designed using Newmark-Hall inelastic 
design response spectra with 5% damping. They c<Hicluded that 
displacement ductility demand on the average is less than the values 
specified by Newmark-Hall method. However, Newmark-Hall method becomes 
less conservative as the values of the specified ductility and viscous 
damping are increased. 

Briseghella et al. (1982) have proposed a relation between M and R^ 
as a function of period based on a study of SDOF systems, for Takeda 



ICO 


model of load-defonnation relationship. They have used artificial time 
histories compatible with the power spectral density of a recorded 
motion d\iring the 1976 Friuli earthquake. 

For EPP load-deformation relationship, Lai and Biggs (1980) have 
investigated the sources of variability of inelastic response spectra, 
due to strong motion duration, ductility, and damping. From time history 
analyses, they concluded that the inelsistic response sp»ectra are not 
significantly dependent on strong ground moticm duration. However, they 
added that this conclusion is valid only vdien the strong ground motion 
with varying durations are compatible with the same prescribed elastic 
response spectrun. They also found that the Newmark-Hall procedure for 
predicting the inelastic response of a 5% damped EET system is 
unconservative, on the other hand for 2% dairping the Newmark-Hall 
approach predicts a conservative inelastic response. Fig. 5.3 shows the 
IKF spectra proposed by them. 

Pal et al . (1987) have developed yield displacement spectra, 
constant strength spectra, constant displacement spectra, reduction 
factor spectra, inelastic acceleration spectra, and inelastic yield 
displacement spectra for three artificial and two actual earthquake 
time-histories using EPP, bilinear, and stiffness degrading models. They 
concluded that the design seismic loads for strong motion earthquake can 
be reduced below those required for elastic analysis by a factor of 
about 1.1 to 15 depending vpcai the time period of the structure and the 
desired displacement ductility. Further, the reduction factors due to 
stiffness degrading model are only slightly higher than those for EPP 

models. This can be seen in Fig. 5.5. 

Miyama et al. (1988) have calculated the inelastic response of 
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SDOF systems for bilinear, and "origin” oriented models of 
load-deformation relationship. The study is based on artificial 
earthquake ground motions of which the target spectnm is elastic 
response spectrum. They have expressed R as functicn of chictility 
factor, natural frequency of systan, and shape of the elastic response 
spectra. 

Riddell et al. (1989) studied non-linear behaviour of EPP SDOF 
system, particularly for short period range, using 53 ground jmoticai 
records, and concluded that reduction factor tends to a value very close 
to unity for very small periods and to a constant value larger than one 
for periods over 0.4 sec. Ihe idealized bilinear relation between and 
T for different values of ductility is shown in Fig. 5.5. 

Using bilinear load-deformation relationship, Miranda (1992) has 
reported a statistical study cxi SDOF systems from 124 grourai motion 
records from different types of site. The- study shows that local site 
conditions significantly affect the value of R^. Fig. 5.7 shows the 
variation of ductility reduction factor with time-period for rock sites 
and for soft soil sites. Notice that while the plot for rock sites is R^ 
versus T, for soft soil the plot is R^ versus T/T^. Here T^ is the 
predominant period of the site defined as the period corresponding to 
the maximum spectral velocity. 

Vidic et al . (1992) state that peaks in the DRF-spectrum correspond 
to peaks in the elastic acceleration spectrum and that, for this reason, 
smooth DRF-spectra are applicable in combination with smooth elastic 
design spectra. Based on the response of an oscillator with Q-hysteretic 
model with 5% danping, they show the influence of ground motions (Fig. 


5.8). 
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Thus, several studies (eg., Riddell et al, 1989, Vidic et al. 1992) 
incltde that the value of 

i) in the short period range depends on T and fi; it increases linearly 
with T, starting with M = 1 and going upto = m, and 
o) in the medium and long period range, depends on #i but is independent 
of T; its value is around /j. 


,.3 METHC® OF MALYSIS 

GRF spectra for all the horizontal component records were obtained 
:or SDOF, EPP systems using Newmark-beta method. The displacemafit 
iuctility value chosen was 2, 3, 5, and 10. The confutations were 
carried out for 41 values of natural periods (0.001, 0.055, 0.01, 0.02, 
0.09, 0.1, 0.2, ..., 3.0). Mean, 

mean-minus-one-standard-deviation, and coeffecient of variation (COV) of 
ductility reduction factors for individual events and for the average of 
six events have been evaluated. 

5.4 RESULTS 

Mean IKF spectra in the six events for displacement ductility of 2, 
3, 5, and 10 are shown in Fig. 5.9. It is seen that, irrespective of the 
value of ductility, for very small periods, DRF is close to unity. 
However, for periods rfto 0.2-0. 3 sec, the spectral values build vp and 
later remain fairly constant upto period of about 2.7 sec. A slight 
increase in DRF value is observed in the period range 2.7 to 3.0 sec. 


5.4.1 Mean DRF spectra 
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The mean 130? spectra for different evmts are compared in Fig. 5.10 
for the four values of ductility. It is seen that for low ductility 
(less than 5) spectral shapes are quite consisteit for different eveits 
with only nordnal variaticais. However, for higher value of ductility the 
variation in mean DRF spectrum tends to increase. Ihis is in line with 
the oteervaticn by Riddell et al. (1989), that "the inelastic response 
of SDCF systems with large inelastic incursions are more affected by the 
characteristics of accelerograams than those with small inelastic 
incursions". 

T^e mean DRF spectra in "Sept" adhieve an musual peak at natural 
period of around 0.2 sec to 0.3 sec. However, beyond 0.4 sec the mean 
DRF spectra of this event are lower than thosefor all the otter events. 
Notice that while the mean elastic spectra of this earthquake (Fig. 
4.16) shows a sharp drop at natural period of around 0.4 sec, the 

elcistic response spectra of this earthquake do not sbsw large 
aimplifications in the period range of 0.2 sec to 0.4 sec. Hence, the WF 
spectra of this evait do not reflect the elastic resp<mse spectrum 
characteristics; this is in contrast with the c4)servations of Vidic et 
al. (1992). The DRF spectra are the highest for evant "April" as 

compared to other events. Interestingly, the anplification factors in 
the acceleration and velocity regicns of elastic spectra for this evaat 
are the lowest. 

5.4.2 OCV and M^ 44 inus-<)ne-Standard-Deviatio^ Spec±raL 

The coeffecient of variation in ordinates of DRF spectra, for 

different ductility, for the six events is shc»m in Fig. 5.11. Ihe 

spectra is shown for period range of 0.1 to 3.0 sec. It is very clear 
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that the OOV is higher for large value of ductility and vise versa. 
Also, irrespective of event, (XN reraains mostly between 0.2 to 0.4, for 
ductility value between 10 to 2. Ihese observations can be clearly seen 
in Pig. 5.12, which shows the average of (XJV spectra from each event. 
Hiis has iapotant inplicaticns for design of important structures where 
the nE?F to be used is more like mean-imnus-one-standard-deviation than 
the mean value. Such a consideration is not being made in the profession 
at present. 

Fig. 5.13 shows the of mean-minus-one-standard-deviation spectra in 
the six events. In general, shape of mean-minus-one-standard-deviation 
spectra is same as that of mean DRF spectra. 


5.4.3 Idealisation of Mean and Itean-^tijiusKJne-Staaidard-Deviaticn Spectra 
Fig. 5.14 and Fig. 5.15 show the mean and 
mean-minus-one-standard-deviation spectra obtained by averaging tlte 
corresponding spectra obtained in the six events. These shapes are 
idealised as bilinear curves and are shown in the sane figure. Similar 
shapes were obtained by Riddell et al. (1989) (Fig. 5.6). The relation 
between E®F and period for idealised shapes of mean and 
mean-minus-one-standard-deviaticn spectra can be expressed as; 





for 0 T < 0.2 sec 
for 0.2 T < 3.0 sec 


(5.6) 


The values of R* for ductility 2, 3, 5, and 10 are givafi in Table 5.1. 
It can be seen from Fig. 5.9, that for longer periods (T > 2.7 sec), for 
high ductility (/^ ^ 5), a slightly higher ductility reductiwi factorthan 



that indicated by idealised shape is affordable. Thus, idealized sha^e 
in this period range is conservative. 

5.5 CCMPARISQN OF RESULTS 

Fig. 5.16 shows the idealized DRF spectra obtained from the data in 
ccrrparison with the DRF spectra proposed by, Newmark and Hall (1973), 
Lai and Biggs (1980) and Riddell et al . (1989) studies. The figure 

clearly indicates that, the variation in different methods exists mostly 
in the period range of 0.2 to about 0.6 sec. The spectra proposed by 
Lai -Biggs is most conservative. The Newmark-Hall spectra niatches quite 
well with the proposed spectra for, periods longer than 0.6 sec. But, 
for shorter periods the Newmark-Hall spectra is quite lower than the 
obtained spectra. The proposed idealised shape of DRf spectra is similar 
to the shape proposed by Riddell et al. (1989). However, there are some 
differences in the period frcam which becomes constant. For higher 
ductility (5 ^ s. 10), Riddell et al . suggests this period as 0.4 sec 

as agairist 0.2 sec proposed from the data. Due to this difference the 
slope of initial linear portion of the two bilinear curves vary. 
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Connparison of ductility reduction factor spectra for 
elasto-plastic and stiffness degrading models (Pal et 
al. 1987). 
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Himalayan earthquakes. 
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CHAPTER VI 

SUMMARY AND CONCLUSIONS 


When formilating design criteria for aseismic design of structures, 
it is necessary to specify the ground motion that the structures must be 
designed to resist. The present knowledge in earthquake engineering and 
related sciences is not enough to know the exact characteristics of the 
future ground motion. However, data available from the past earthquakes 
at the site could be effectively used to estimate the earthquake ground 
motion in future earthquake. 

In the present study, an attempt has been made to study the data 
available from the six Himalayan earthquakes, which occurred during 
1986-1991. A data set consisting of 198 horizontal component and 99 
vertical component records is used to study the ground motion 
characteristics such as peak parameters, duration of strong ground 
motion, shape of response spectrum, and ductility reduction factor. The 
results obtained were compared with those available from similar studies 
on earthquakes elsewhere, and with currently adopted design practices. 

Following are the salient conclusions regarding characteristics of 
ground motions in the Himalayas from the six earthquakes recorded in the 
region: 

- The A/V ratio for the ground motion in the Himalayan region is 
unusually high. The mean A/V is more like 20.0 sec for rock sites, 
this is quite high when compared to the value of 12-15 sec ^ reported 
for earthquakes elsewhere. The mean-minus-one-standard-deviation A/V 
value is obtained as 12 sec 

- The variation in A/V ratio with the distance from shock-source is 



sonewhat insignificant for the Himalayan earthquake ground imtions. 
Thus, high A/V value is obtained even at large distances. This 
suggests that attenuation characteristics of peak ground motions in 
this region is different from those reported in literature, 

- The mean and mean-plus-one-standard-deviation value of AD/V^ are 

obtained as 4.6 and 8.0, respectively. No clear trend in variaticn of 
2 

AD/V with distance is observed. 

- The mean and mean-plus-one-standard-deviation value is datained 

as 0.51 and 0.70, respectively. A conservative value of 0.67 or 0.75 
is suggested as against the currently adopted value of 0.5 in Indian 
code. 

- A/V is found to be good indicator of frequency conteat of ground 
motion and hence the relative shape of spectra. 

- Duration of strong shaking, as evaluated by Trifunac and Brady (1975) 
definition, shows no definite relationship with A/V or with distance 
from shock-source, though such relationships are reported in 
literature. 

- Significant variation in elastic spectral shapes in intermediate 
period range are observed depending on the magnitude of earthquake. 
Mean elastic spectral shapes for three ranges of earthquake 
magnitude; (1) small magnitude (M < 5.5), (2) moderate magnitude (5.5 
i M i 6.5), and (3) large magnitude (M > 6.5) are obtained. This could 
be helpful in deciding the shape of design spectra. 

- Spectral shapes do not vary much with distance for moderate magnitude 

4 . 

earthquakes within moderate epi central distances. 

- Both, the ratios of peak parameters and the elastic response spectra 
indicate that with distance the attenuation of waves of different 
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frequencies is not so different as elsewhere. Hence, for projects in 
Himalayas, the empirical relationship for peak parameters with 
magnitude, epicentral distance, etc. , of elsewhere are not applicable. 
The design spectral shape provided in UBC (1991) is quite conservative 
if applied to Himalayan earthquake ground motion as it ccnpletely 
envelops the obtained mean shape. The shape of design spectra 
currently adopted in 13:1893-1984 is quite different from the mean 
shape obtained frctn data. The code spectra underestimates the spectral 
ordinates in period range of 0.1 to 0.35 sec and overestimates in the 
period range greater than 0.35 sec. Revision of currently adopted 
shape in Indian code is urgently required. In order to bring the mean 
spectral shape of Indian code at par with the shape obtained frcm the 
data, the acceleration amplification in the period range of 0.1 to 
0,35 sec may be increased, while the presently used shape for longer 
periods is conservative and can be kept as it is. 

Shapes of AERB "standard spectra” for rock sites and that of 
mean-plus-one-standard-deviation spectra from the data are in good 
agreannent. This indicates that spectral shapes suggested by AS® are 
quite satisfactory for the structures which are usually designed using 
mean-plus-one-stahdard-deviation spectra . 

The acceleration amplification region mostly ranges between natural 
periods of 0.12-0. 26 sec, while velocity region ranges between periods 
0.26 to 2.6 sec. Mean acceleration amplification factor for the data 
is higher than that for data from earthquakes elsewhere. Qn the other 
hand, mean amplification factor in velocity region is lower. However, 
mean-plus-one-standard-deviation acceleration amplification from 
Indian data is comparable with those obtained in other studies. 
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Straight line correlation between spectrum intensity and peak ground 
velocity are established. It is found that the slope of straight line 
correlating spectrum intensity with peak ground velocity increases 
with the increase in earthquake magnitude, 

- A relationship between ductility reduction factor and time period of 
elasto-plastic single degree of freedom structure is proposed. Except 
for some difference for periods less than 0.4 sec the proposed spectra 
is quite similar to the one proposed by Riddell et al. (1989). 

Based on the preset work, following area for research in future are 
identified. 

- A/V is an excellent indicator of frequency content. Hence, design 
recorrmendation can be based on A/V ratio. This approach is adopted by 
National Building Code of Canada (NBCC 1990). Hence, study of A/V 
ratio to work out design recommendations can be carried out. 

- Anplification factors in displacement region could not be evaluated 
because of the unexpected shape of spectra in this region. A detailed 
investigation on this is needed. 

- The present work on ductility reduction factor is limited to 
elasto-plastic ■ single degree of freedom systems, this work could be 
extended for other load-deformation models and for multi degree of 
freedom systems. 

- Now better strong motion instrumentation is available, to make full 
use of it the ground motion parameters like magnitude, location of 
epicentre, etc., are needed with better precision. This will be very 
helpful in developing the attenuation characteristics of the peak 
ground motion parameters. 

- In India, detailed studies on analysis of available strong moticai data 
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are only few, so many other parameters and characteristics could be 
studied so that cme could reach to the current state-of-art in 
California or Japan. 
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APPENDIX A 


Table A-l For April 26, 1986 earthpuake recorded in Kanpra Arrav. 




EPICENTRAL 





DURATION 

'spectruR^ 

SR, 

LOCATION 

DISTANCE 

dirn! 

PEAK 

PEAK 

PEAK 

TRIFUNAC 

INTENSITY 

NO, 





GROUND 

GROUND 

GROUND 

BRADY 

FOR 5% 






ACCL. 

VEL. 

DI5PL. 

(1975) 

DAMPING 




(Km) 


(q ) 

( mm/s) 

(mm ) 

(sec ) 

(mm ) 

1 

BHANDLKHAS 

A1 

= 14.8 

L 

0.145 

83.13 

20.08 

2,82 

174.92 



A2 

= 22.3 

T 

0.125 

90.75 

5.85 

3.16 

147.23 



A3 

= 26.3 

V 

0.020 

19.03 

6.33 

9.04 

61.93 



A4 

= 24.9 







1 

BAROH 

A1 

» 10.1 

L 

0.058 

36.80 

6.55 

4.34 

09.03 



A2 

= 13.1 

T 

0.057 

27.00 

4.58 

6.74 

101.74 



A3 

= 22.2 

V 

0.020 

15.07 

3.64 

10.50 

57.85 



A4 

= 23.2 







3 

BHAWARNA 

A1 

= 14.7 

L 

0.036 

11.70 

2.73 

9.04 

58.39 



A2 

= 19.1 

T 

0.035 

19.21 

4.20 

8,09 

75.72 



A3 

= 26.3 

V 

0.036 

9.34 

21.13 

5.98 

36.85 



A4 

= 25.5 







4 

DHARANSHALA 

A1 

= 9.9 

L 

0.176 

72.97 

7.76 

2.60 

188.37 



A2 

= 11.3 

T 

0.186 

94,90 

24.79 

3.94 

379.46 



A3 

= 5.0 

V 

0.0B3 

27,39 

4.19 

4.34 

99.61 



A4 

= 2.8 







5 

JAWALI 

A1 

= 36.5 

L 

0.015 

19,20 

6.90 

15.92 

70.62 



A2 

= 29.0 

T 

0.017 

15.18 

9.56 

13.99 

58.34 



A3 

= 26.4 

V 

0.011 

14.39 

7.53 

16.14 

71.01 



A4 

= 28.6 







6 

KANGRA 

A1 

= 13.5 

L 

0.148 

50,57 

5.23 

4.84 

143.84 


A2 

a 5.2 

1 T 

0. n 1 

95 . 75 

^.6? 

3.74 

740 . 7/t 



A3 

= 11.0 

V 

0.072 

31.93 

5.37 

5.40 

104.75 



A4 

= 12.5 

j 





228 . 38 

7 

NABROTABAGWAN 

A1 

« 3.6 

! L 

0.149 

94.17 

13.09 

4.20 


A2 

= 6.4 

T 

0.080 

25.40 

6.14 

5.14 

78.55 



A3 

= 12.9 

V 

0.050 

18.19 

6.51 

14.82 

89.39 



A4 

= 12.5 






193.59 

309.97 

8 

SHAHPUR 

A1 

A2 

= 20.4 
= 15.8 

L 

T 

0.200 

0.250 

59.21 
147.80 1 

7.26 

10.85 

2.70 ! 

2.00 i 



A3 

= 6.8 

V 

0 . 065 j 

28.04 

5.25 

6.40 i 

110.13 



A4 

= 10.3 




j 



9 

SIHUNTA 

A1 

A2 

= 34.3 
= 30.6 

1 1 

1 L 
‘ T 

0.050 

0.036 

28.72 

33.88 

4.71 

3.67 

6.20 ' 
U m 66 

96.32 

105.27 



' A3 

= 22.4 

V 

0.039 

27.73 

4.74 

9.34 : 

78.79 



A4 

= 23.4 





i 


L ; 

Lonq i tudinal , 

L, — 

T : 

Transverse, 

V : Vertical 
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Tablff A-2 For Bfpttmbtr 10, 1906 ©trthquak# recorded in Shi 1 long array. 




“TFTCrNTRSr 





[DURAflON 

IspecYrum 

SR« 

LOCATION 

DISTANCE 

DIRN. 

PEAK 

PEAK 

PEAK 

TRIFUNAC 

INTENSITY 

NO, 





GROUND 

GROUND 

GROUND 

Zi BRADY 

FOR 5% 






ACCL. 

VEL. 

DISPL. 

(1975) 

DAMPING 




(Km) 


(g) 

(mm/s ) 

(mm ) 

(sec ) 

(mm) 

1 

BAITHLANGSO 

Al 

» 139.4 

L 

0.045 

20.58 

3.55 

10.38 

60.31 



A2 

= 102.1 

T 

0.042 

12.07 

2.73 

10.78 

55.40 



A3 

= 181.2 

V 

0.025 

10.91 

4.74 

11.86 

39.24 

2 

DAUKI 

Al 

a 42.7 

L 

0-089 

32.40 

3.00 

3.92 

64.50 



A2 

a 32.5 

T 

0.090 

37.33 

2.56 

5.30 

67.49 



A3 

a 48.1 

V 

0.032 

14.14 

5.92 

14.26 

60.31 

3 

KHLIEHRIAT 

Al 

a 01.6 

L 

0.031 

11.30 

2.00 

10.56 

35.00 



A2 

a 56.9 

T 

0.046 

18.41 

3.86 

9.60 

49.81 



A3 

a 54.4 

V 

0,001 

7.71 

1.95 

11.08 

37.01 

4 

NONGKHAW 

Al 

a 65.4 

L 

0.055 

33.09 

9.32 

17.24 

103.41 



A2 

a 35.8 

T 

0.093 

48.08 

5.47 

11.76 

120.43 



A3 

a 28.3 

V 

0.034 

14.11 

6.83 

23.50 

64.77 

5 

NONGPOH 

Al 

a 93.5 

L 

0.054 

21.13 

4.26 

8.32 

39.04 



A2 

a 56.4 

T 

0.056 

11.32 

4.59 

12.86 

35.94 



A3 

a 33.8 

V 

0.033 

11.25 

2.84 

11.19 

27.90 

6 

NONGSTOIN 

Al 

a 55.8 

L 

0.019 

10.04 

4.90 

7.76 

32.69 



A2 

54.5 

T 

0.014 

6.4e 

i .66 

7.!5'J 1 




A3 

a 64.3 

V 

0.008 

5.95 

1.99 

7.82 

22.68 

7 * 

PAN I MUR 

Al 

a 134.9 

1 

' L 

0,040 

i 9.37 

2.02 

8.94 

27.08 



A2 

a 103.8 

T 

! 0.048 

21.15 

' 2.65 

5,78 

39,76 



A3 

a 90,2 

V 

0.023 

5,61 

1.40 

9.42 

27,87 

8 

PYNURSLA 

Al 

a 37.4 

L 

0,093 

26.75 

3.24 

5.14 

83.82 



A2 

a 15.0 

T 

0.076 

20.37 

4.55 

8.68 

62.32 



A3 

a 33.6 

V 

0.030 

9.28 

3.80 

14.78 

35.17 

9 

SAITS AM A 

Al 

= 104.4 

L 

0,110 

36.69 

4.90 

6.12 

74,97 


A2 

a 68.3 

T 

0.140 

58.73 

4.92 

2.44 

106.37 



A3 

a 50.1 

V 

0.062 

20.06 

3,22 

7.32 

54,22 

10 

UMMUL0N6 

Al 

a 70.9 

L 

0.110 

26.53 

2.38 

5.82 

46.13 

A2 

a 36.7 

T 

0,063 

12.10 

1.56 

9.56 

34.91 



A3 

a 27.4 

V 

0,050 

14.25 

3.21 

11.42 

47.99 

11 

UMRONGSO 

Al 

=111.1 

L 

0,027 

9.68 

1.30 

8.26 

32.03 

A2 

a 02.4 

T 

0.032 

11.84 

4.41 

9.82 

35.96 



A3 

a 72.5 

V 

0.014 

7.81 

2.26 

10,66 

32.72 

12 

UMSNING 

Al 

A2 

a 76.4 
a 38.3 

L 

T 

0,100 

0.076 

28,72 

25.67 

3.42 

3.93 

4.26 

8.78 

57.56 

64.84 



A3 

a 15.0 

V 

0.050 

10.88 

3.75 

14,50 

41.64 


# L t Longitudinal, T i Transv,rs». V t VBrtical 


Table A 3 For May IS^j 19S7 earthquake recorded in Shillong array* 




EPICENTRAL 





DURATION 

SPECTRUM 

SR. 

LOCATION 

DISTANCE 

DIRNt 

PEAK 

PEAK 

PEAK 

TRIFUNAC 

INTENSITY 

NO. 






GROUND 

GROUND 

GROUND 

Zi BRADY 

FOR 5% 







ACCL. 

VEL. 

DISPL . 

(1975) 

DAMPING 




(Km) 


(g) 

(mm/s ) 

( mm ) 

(sec ) 

{ mm ) 

1 

BAITHLANGSO 

A1 


121.8 

L 

0.034 

15.53 

4.15 

16.84 

55.89 



A2 

= 

177.5 

T 

0.027 

21 .56 

8.11 

18,50 

70.02 



A3 

ss 

112.8 

V 

0.020 

12.85 

3.80 

18.14 

59.79 

2 

BAMUN6A0 

A1 

sr 

82.1 

L 

0.020 

12.22 

4.24 

22.58 

44.31 



A2 

S! 

137.7 

T 

0.020 

22.00 

6.78 

23.48 

57.87 



A3 


73.8 

V 

0.019 

13.62 

3.32 

25.18 

65.10 

3 

BERLONGFER 

A1 


54.0 

L 

0.072 

35.95 

5.51 

26.52 

110.19 



A2 

= 

109.5 

T 

0.090 

36.20 

8.95 

22.78 

128.03 

1 


1 

s 

46 .-1? 

V 


19.24 

7.47 

30.62 

72.52 

4 

BOKAJAN 

A1 

2= 

56.4 

L 

0.030 

25.96 

9.34 

18.88 

75.47 



A2 

= 

91.5 

T 

0.066 

34.28 

10.31 

12.04 

107.65 



A3 


58.4 

V 

0.020 

13.20 

5.25 

21.92 

44.61 

5 

DIPHU 

A1 

ss 

53.7 

L 

0.086 

27.69 

5.80 

27.06 

94.39 



A2 

S 

104.7 

T 

0.073 

23.92 

11.54 

27. 18 

90.79 



A3 

= 

49.6 

V 

0.054 

15.31 

6.65 

26.68 

78.30 


fiUNJUNfi 


2= 

73.5 

L 

0.040 

23.04 

5.55 

12.46 

84.34 



A2 

S 

120.0 

T 

0.049 

1 28.56 

5.51 

11,60 

i 90.31 



A3 

= 

64.0 

V 

I 

0.019 

10.83 

3.38 

1 12.08 

I 

! 53.12 

1 

1 

7 

HAFLONG 

A1 

s 

77.9 

t 

L 

0.055 

^ 36.93 

8.35 

1 

6.62 

120.59 


A2 


119.0 

1 T 

0.035 

20.61 

5.38 

9.82 

107.89 


1 

A3 

= 

69.3 

V 

1 0.015 

j 

1 12.15 

i 

3.06 

11.06 

62.22 

8 

1 

HAJADISA 

|A1 

=: 

42.7 

L 

I 

0.078 

i 38.07 

7,50 

12.16 

111.86 

A2 

- 

90,5 

T 

0.086 

28.07 

5.00 

11.24 

99 . 50 



1 ^ 

1 A3 

= 

33.6 

V 

0.028 

1 

15.34 

3,49 

13.40 

74.82 

9 

HATIKHAL I 

i 

1 

1 A1 


62.3 

i 

L 

1 

0.031 

14.55 

5.55 

27.06 

60.68 


i 

: A2 


117.5 

T 

i 0.038 

17.71 

3.48 

27.26 ; 

70 . 48 



A3 

= 

52.7 

V 

0.027 

18.58 

3.64 

28.86 

74.64 

10 

LAISONG 

1 

Al 

A2 


52.2 

90.4 

L 

T 

0.042 
; 0.061 

32.25 

23.69 

7,84 

9.42 

11.28 

9.68 

103.30 
109. 10 



A3 

= 

45.0 

V 

0.019 

13.54 

6.20 

13.68 

67.30 

11 

NONGPOH 

Al 

A2 

A3 

= 

188.7 

243.4 

L 

T 

0.017 

0.047 

13.12 

12.58 

4.88 

5.93 

16.24 

16.46 

35.71 

51,92 



“ 

179.0 

V 

0.013 

1 9.55! 

1 

3. 13 

16.62 

54.64 

2 

PANIMUR 

Al 

A2 

A3 

a 

91.6 

146,0 

L 

T 

0.040 
i 0.047 

1 

18.00 

1 18.11, 

4. 10 
2.99 

9. 12 
9.92 

53.32 

47.74 



= 

81.7 

V 

; 0.017 

1 16.791 

5.78 

9.52 

37 . 23 

. con td . 







136 


SR. 
NO. i 

LOCATION 

EPi CENTRAL 
DISTANCE 

(Km ) 

dirnT 

PEAK ! 
GROUND 
ACCL. 

(Q) 

i 

PEAK 

GROUND 

VEL. 

(mm/s > 

PEAK 

GROUND 

DISPL. 

(mm ) 

''DURAT'IW ' 
TRIFUNAC 
& BRADY : 

(1975) 

(sec) 

-SFEcrmfr' 

INTENSITY 
FOR 5% 

DAMPING 

( mm ) 

13 

SAITSAMA 

A1 

- 

1-33 . 9 

L 

0.037 

13.92 

3.36 

20 . 36 

39. 16 



A2 

= 

188. 1 

T 

0.050 

22.84 

8.91 

20 . 38 1 

80. 19 



A3 


124. 1 

V 

0.019 

9.29 

4.00 

21.14 ; 

45 . 12 

14 

UhRONGSO 

A1 


106.3 

L 

0.021 

13.79 

2.76 

10,94 

44.37 


[ 

A2 

= 

160.2 

T 

0.026 

14.50 

5.69 

9.90 

51 .81 


L- 

|A3 

= 

98.3 

V 

1 0.016 

9.13 

2.91 

10.04 

41.11 


L : Longitudinal, T : Transverse, V ; Vertical 
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SR- LOCATION 
NO. 


1 BAISAO 


BAITHLANGSO 


5 BAMUNGAO 


DAUKI 


GUN JUNG 


HAFL0N6 


HATIKHALI 


KATAKHAL 


KHLIEHRIAT 


NAWFHLANG 


N0N6KHLAW 


NONGPOH 




* 




DURATION 

SPECTRUM 

UiblWNCE 

DIRN. 

PEAK 

PEAK 

PEAK 

TRIFUNAC 

INTENSIT 




GROUND 

GROUND 

GROUND 

h BRADY 

FOR 5% 




ACCL, 

VEL, 

DISPL . 

(1975) 

DAMPING 


(Km) 


(g) 

(mm/s ) 

(mm ) 

(sec ) 

(mm ) 

A1 

= 190.3 

L 

0.022 

14.14 

3.90 

11,28 

56, 19 

A2 

= 159.0 

T 

0.024 

20.64 

6.90 

10.72 

115.95 

A3 

= 150.4 

V 

0.009 

1 1 .27 

3.26 

10,58 

73.46 

A1 

= 193.2 

L 

0,031 

13.65 

3.81 

10,06 

59,99 

A2 

= 184.0 

T 

0.022 

12.26 

3.33 

9,94 

70.49 

A3 

= 135.5 

V 

0,014 

7.44 

1.76 

10.46 

29.40 

A1 

= 223.1 

L 

0.016 

5.08 

1,30 

7.24 

27.01 

A2 

= 203.9 

T 

0.013 

10.35 

2.30 

7.72 

33.88 

A3 

= 175.5 

V 

0.010 

6.23 

2.33 

7.20 

26.98 

A1 

= 103.7 

L 

0.027 

18.01 

2.50 

7.50 

55.74 

A2 

= 79.0 

T 

0.038 

18.37 

4.41 

6.78 

57.38 

A3 

= 75.5 

V 

0.010 

8.57 

3.41 

8.62 

33.73 

A1 

= 202.8 

L 

0.036 

29,05 

6.26 

8.68 

103.77 

A2 

= 166.7 

T 

0.037 

24.44 

3.87 

7.96 

74.97 

A3 

= 166.5 

V 

0.020 

12.87 

3.30 

10.78 

44,54 

A1 

= 202.5 

L 

0.034 

19.42 

3. 13 

8. 16 

67.35 

A2 

= 162.3 

T 

1 0.027 

14.02 

! 3.34 

8.04 

71 ,51 

A3 

= 170.8 

V 

0.008 

7.22 

1 .99 

9,04 

22.00 

A1 

= 222.3 

L 

0.024 

8.78 

3.18 

8.98 

33.11 

A2 

= 195.3 

T 

0.025 

9,87 

3,84 

10. 18 

33.50 

A3 

= 175.8 

V 

0.021 

7,29 

1 .98 

10.34 

28.00 

A1 

= 164,6 

L 

0,009 

10,64 

3.04 

8,82 

64.70 

A2 

= 114,9 

T 

0.008 

7,75 

1 .60 

9,28 

32,11 

A3 

= 149.4 

V 

0.010 

9,88 

4.00 

9.30 

44.81 

A1 

= 142,2 

L 

0.080 

47.46 

8.45 

3,90 

133.52 

A2 

= 116.9 

T 

0.066 

27.87 

3,44 

6. 12 

85,16 

A3 

= 103.0 

V 

0.029 

15.91 

2,94 

9.30 

49,40 

A1 

= 91.7 

L 

0.081 

43.85 

7.97 

19.90 

98. 18 

A2 

= 93.5 

T 

0,059 

24,89 

10.14 

20.28 

100.55 

A3 

= 42.0 

V 

0,036 

12.33 

4.86 

23.74 

42,19 

A1 

= 99.3 

L 

0,110 

42.94 

9.58 

24.20 

89.75 

A2 

= 116.5 

T 

0.114 

52.63 

11.41 

13,52 

156.09 

A3 

= 35.0 

V 

0.104 

35.10 

11.23 

16.58 

84.15 

A1 

132.9 

L 

0.027 

15.89 

4.64 

12.28 

46.50 

A2 

144.7 

T 

0.086 

26,45 

3.75 

8.40 

74.80 

A3 

= 68.4 

V 

0.038 

10.64 

4.02 

12.68 

40. 16 


con td 
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SR. 

NO. 

LOCATION 

epicentral 

DISTANCE 

(Km) 

dirnT 

PEAK 

GROUND 

ACCL. 

(g) 

PEAK 

GROUND 

VEL. 

(mm/s ) 

PEAK 

GROUND 

DISPL. 

(mm ) 

DURATION 
TRIFUNAC 
Si BRADY 
C1975) 

Csec ) 

SPECTRUM 
INTENSITY 
FOR 57. 
DAMPING 

(mm) 

13 

PYNURSLA 

A1 

= 

96.1 

■■ 

0.050 

21.24 

4.99 

21.96 

72.09 



A2 


83.1 


0.031 

11.95 

4.93 

22.54 

51.90 



A3 


59.4 

H 

0.015 

9.43 

3.63 

28.44 

47.79 

14 

SAITSAMA 

A1 


159.9 

■ 

0.066 

22.38 

4.50 

9.28 

71.83 



A2 

= 

148-6 

T 

0.058 

32.43 

4.37 

9.66 

75.86 



A3 

= 

106. 1 

V 

0.032 

9.90 

2.37 

12.52 

48.50 

15 

SHILLONG 

A1 

= 

109.1 

L 

0.047 

16.06 

5.16 

8,82 

57.06 



A2 


108.4 

T 

0.036 

15.39 

7.86 

9.00 

64.78 


1 

A3 


55.3 

V 

0.013 

14.42 

[ 

5.19 

1 

; 10.64 

57 . 49 

16 

UmULONG 

A1 

= 

127.9 

L 

0.056 

22.48 

^ 5.69 

19.82 

66 . 66 



A2 

s 

115.3 

T 

0.055 

21.42 

6.35 

19.18 

81.74 



A3 

= 

80.0 

V 

0.024 

10.89 

4.37 

22.30 

40,30 

17 

UMRONGSO 

A1 

ss 

171 . 1 

L 

0.046 

31.43 

3.06 

7.80 

81 .06 



A2 

= 

146.7 

T 

0.037 

26.58 

6.67 

9.32 

95.65 



A3 

r 

126.9 

V 

0.022 

12.26 

2.38 

11.20 

47.64 

8 

UfiSN I NG 

A1 

= 

205.5 

L 

0.040 

16.62 

3. 14 

18. 19 

51.45 



A2 

- 

184.0 

T 

0.061 

36.93 

5.03 

14.54 

78.16 



A3 


155-5 

V 

0.018 

11.95 

3.55 

20.84 

53.25 


L : Longitudinal, T : Transverse, V : Vertical 
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fable A 5 For August 6, 1988 earthquake recorded in Shillong array 



LOCATION 

EP I CENTRAL 





DURATION 

SPECTRUM 

SR. 

DISTANCE 

DIRN. 

PEAK 

PEAK 

PEAK 

TRIFUNAC 

INTENSITY 

NO. 






GROUND 

GROUND 

GROUND 

BRADY 

FOR 5% 







ACCL. 

VEL. 

DISPL. 

(1975) 

DAMPING 




(Km> 


(g) 

<mm/s ) 

<mm ) 

( sec ) 

(mm ) 

1 

BAIGAO 

A1 

=: 

227.5 

L 

0.221 

65.68 

10. 19 

31,919 

169.00 



A2 


226.7 

T 

0. 144 

63.70 

9.80 

33.899 

182,74 



A3 

= 

165.0 

V 

0.053 

20.11 

5.85 

38.180 

83.00 

2 

BAITHLANGSO 

A1 


268.2 

L 

0. 154 

75.52 

12.36 

28.079 

256.94 



A2 

= 

266.2 

T 

0. 165 

117.05 

15.45 

28.839 

420 . 49 



A3 


200.3 

V 

0.082 

36.12 

6.98 

46.099 

136.21 

3 

BAMUNGAO 

A1 

s: 

227.6 

L 

0.093 

64.03 

8.96 

24.319 

203.34 



A2 

= 

225.5 

T 

0.069 

49.80 

10.23 

22 . 979 

200.59 



A3 


159.5 

V 

0.065 

25.06 

4.23 

23.359 

85.14 


BERLONGFOR 

A1 

=: 

199,5 

L 

0.300 

217.24 

33.40 

21.960 

898 . 1 1 



A2 

=: 

197.5 

T 

0.344 

228.19 

36.30 

20 . 000 

806.13 

■ 


A3 

= 

131.6 

V 

0.174 

90.26 

13.18 

41.599 

247.96 

5 

BOKAJAN 

A1 


167.2 

L 

0. 151 

86.75 

19.91 

20 . 699 

345 . 83 



A2 

=: 

163.8 

T 

0.224 

121.36 

20.41 

20.499 

369.27 



A3 

= 

99.5 

V 

0. 148 

32.49 

9, 14 

30.819 

141.69 

> 

CHERAPUNJI 

A1 

s; 

337.2 

L 

0.052 

21.47 

3,57 

15.839 

92.27 



A2 

- 

336.9 

( T 

i 0,055 

26.78 

6 . 22 

15,219 

111,90 



A3 


277.0 

: V 

0.024 

20.45 

3.81 

17.780 

i 65.50 

i 


DAUKI 

A1 


308.5 

L 

0.110 

46.87 

7.56 

1 

15.779 

i 152-82 



^A2 


308.4 

T 

0.073 

44.53 

5.78 

' 14.899 

1 150.53 



A3 

= 

249.5 

V 

0.031 

20.74 

6.52 

: 18,139 

73.07 


^DIPHU 

A1 

sr 

190.1 

L 

0.283 

181 .48 

23.25 

19.979 

i 636.62 



A2 

= 

187.5 

T 

0.338 

205.55 

22.74 

16,779 

' 561,12 



A3 

= 

121.4 

V 

0, 180 

56.12 

9,00 

34.659 

196.14 


DUlOO 

A1 


232.6 

L 

0.064 

58.89 

14.40 

18,739 

: 274.51 


A2 

= 

233.1 

T 

0.062 

53.49 

12.61 

20,099 

238 . 49 



A3 

= 

180.0 

V 

0.039 

33.83 

7.63 

28.079 

• 156.03 

) 

GUNJUNG 

A1 

_ 

211.2 

L 

0.092 

44.19 

8.74 

28. 199 

176,73 


A2 


210.7 

T 

0.133 

52.05 

8,58 

25 , 239 

202,44 



A3 

= 

150.4 

V 

0,062 

25.56 

6.74 

41 .279 

88.27 


HAJADISA 

A1 


182.5 

L 

0.092 

42.73 

9.72 

34, 180 

175.12 


A2 


181.6 

T 

0.098 

46.43 

8.77 

32.180 

174.97 



A3 


120.1 

V 

0.046 

22.28 

5.98 

43.599 

89,18 

1 

HARENGJAO 

A1 

A2 

A3 


225.0 

225.1 

L 

T 

0.065 

0.078 

40.99 

44.39 

7.57 

8,92 

20,839 

19,459 

167.09 

161.06 




168.4 

V 

0.032 

22. 15 

6.31 

21,260 

69.89 
, . con td . 
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SR« 

NO. 

LOCATION 

EPICENTRAL 

DISTANCE 

(Km) 

dirnT 

PEAK 

BROUND 

ACCL. 

( 9 ) 

PEAK 

GROUND 

VEL. 

(mm/s ) 

PEAK 

GROUND 

DISPL, 

(mm) 

DURATION 
TRIFUNAC 
h BRADY 
(1975) 

(sec ) 

SPECTRUH 
INTENSITY 
FOR 5% 
DAMPING 

(mm ) 

13 

HOJAI 

A1 

= 245.8 

L 

0.108 

51.88 

14.06 

22.439 

203.24 



A2 

= 243.6 

T 

0.134 

65.61 

11.61 

21.939 

260 . 09 



A3 

= 177.5 

V 

0.061 

25.25 

7.83 

41.900 

121,04 

14 

JALLALPUR 

A1 

= 265.2 

L 

0.029 

30.40 

4.81 

13.559 

130,38 



A2 

= 265.5 

T 

0.023 

20.46 

5.18 

13.879 

82.49 



A3 

= 209.8 

V 

0.015 

17.07 

4.01 

13.419 

57.62 

15 

JHIRI 6 HAT 

A1 

202.8 

L 

0.098 

88.30 

11.15 

18.500 

364.48 



A2 

= 203.8 

T 

0.089 

70,35 

10.05 

20.779 

296.53 



A3 

= 155.3 

Y 

0.031 

30.38 

8.21 

25.700 

148, 15 

16 

KALAIN 

A1 

= 253.8 

L 

0.057 

77,76 

15.87 

17.979 

328.96 



A2 

= 254.2 

T 

0.051 

51.09 

15,71 

19.799 

223.70 



A3 

= 199.3 

V 

0.029 

43.21 

10.89 

22.620 

147,82 

[7 

katakhal 

A1 

= 248.7 

L 

0.067 

110.74 

23.34 

18.360 

506,98 


1 

A2 

= 249.5 

I T 

0.059 

I 93.56 

: 25.26 

19.339 

366.33 



A3 

= 197.8 

1 ^ 

0.018 

1 33.72 

10.11 

21 . 199 

1 

133.94 

8 

KHLIEHRIAT 

A1 

= 274.8 

L 

0.070 

29.21 

4.55 

26.780 

105.75 



A2 

= 274.2 

T 

0.071 

31,74 

5.42 

31 .659 

112.44 



A3 

= 213.4 

V 

0.035 

16.1 1 

3,44 

46.219 

47,66 

9 

KOOHBER 

Ai 

= 210.7 

L 

0.049 

39,04 

9.86 

1 6 . 899 

184.96 



A2 

= 211.2 

T 

0.037 

34.88 

7. 14 

18,920 

169,30 



A3 

= 158.2 

V 

0.027 

23.23 

5.57 

19.920 

98.26 

0 

LOHARGHAT 

Al 

376.1 

L 

0.058 

36.58 

8.58 

30,879 

205 . 47 



A2 

= 374.1 

T 

0-054 

48,61 

9.91 

28.459 

255. 16 



A3 

= 303.4 

V 

0,022 

20.42 

5.06 

27.360 

94.73 

1 

HAWKYRWAT 

Ai 

= 364,6 

L 

0.046 

25.18 

3, 10 

17.799 

70,94 



A2 

= 364.1 

T 

0.046 

21 ,98 

4.80 

18.920 

78.82 



A3 

= 303.4 

V 

0.032 

17,51 

4.69 

18.360 

66,87 

7 

HAWPHLANG 

Al 

= 335.6 

L 

0.115 

39,17 

6.78 

25.359 

102.96 



A2 

= 335.0 

T 

0. 107 

38,99 

5.34 

29.659 

83.51 



A3 

= 273.4 

V 

0.036 

12,52 

3,85 

30,019 

45.56 


. . contd . 
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SR. 

NO. 

LOCATION 

EPI CENTRAL 
DISTANCE 

(Km ) 

DIRN* 

PEAK 

GROUND 

ACCL. 

(Q> 

PEAK 

GROUND 

VEL. 

(mm/s ) 

PEAK 

GROUND 

DISPL. 

(mm ) 


SPECTRUM 
INTENSITY 
FOR 5% 
DAMPING 

(mm ) 

23 

MAWSYNRAM 

A1 

= 352.3 

L 

0.085 

33.67 

8.29 

17.399 

136.06 



A2 

= 352.0 

T 

0.064 

28.61 

5.47 

18. 199 

124.34 



A3 

= 292.0 

V 

0.036 

21.45 

5.89 

17.760 

90.86 

24 

NONGKHLAW 

A1 

= 352.8 

L 

0.138 

64.50 

6,02 

30.879 

149.06 



A2 

= 351.8 

T 

0.146 

55.66 

11,07 

28.100 

182.58 



A3 

= 288,4 

V 

0.083 

34.17 

9.46 

31 .399 

97.47 

25 

NONGSTOIN 

A1 

= 385,9 

L 

0.054 

24.09 

9.54 

28.399 

83.96 



A2 

= 385.3 

T 

0.052 

31.04 

11.16 

27.779 

99,32 



A3 

= 323,5 

V 

0.042 

14.80 

6.14 

29,540 

65.43 

26 

PAN I MUR 

A1 

= 238,1 

L 

0.169 

55. 19 

5.29 

31.379 

132.72 



A2 

= 236,7 

T 

0. 125 

44.38 

6.06 

41.819 j 

110.47 



A3 

= 172.3 

V 

0.073 

20.85 

6. 14 : 

48.159 

59.68 

27 1 

PYNURSLA 

A1 

= 320,8 

L 

0.049 

42.13 

5.92 

22 . 439 

j 138.09 



A2 

= 320.4 

T 

0,051 

32.24 

7,22 

24.200 

: 112.20 



A3 

1 

= 260,2 

V 

0.036 

31,55 

j 6.35 

23.159 

116.57 

Z8 

SAITSAMA 

1 

Ai 

= 280.3 

L 

0.211 

' 87.31 

7.83 

33,799 

218.48 


i 

A2 

= 279.0 

T 

0.233 

103.63 

13.89 

32.519 

261 .50 



A3 

= 214,6 

V 

0.099 

33.77 

4.82 

46.840 

95,51 

29 

SHILLONG 

Al 

= 323.9 

L 

0.075 

20.44 

3.24 

24,319 

65.13 



A2 

= 323,1 

T 

0.057 

17.55 

5.86 

26. 139 

58.26 



A3 

= 260,5 

V 

0.036 

9.75 

3.31 

26,000 

42.02 

50 

SILCHAR 

Al 

= 233,0 

L 

0,064 

71,52 

21 . 10 

22.579 

277.62 



A2 

= 233.0 

T 

0,091 

98.76 

21 , 13 

18.079 

426,91 



A3 

= 188-1 

V 

0.025 

29.98 

8.62 

41.720 

98.40 

n 

UMMULONG 

Al 

= 298.4 

L 

0.100 

30.76 

7,86 

30,639 

87.60 



A2 

= 297.6 

T 

0.153 

43.95 

12.80 

25 . 780 

128.93 



A3 

= 235.3 

V 

0.062 

27,33 

11.28 

34.859 

88.15 

:2 

UMRONGSO 

Al 

= 252.3 

L 

0,780 

51.66 

6.46 

31.359 

197,65 



A2 

= 251.4 

T 

0.081 

39,98 

7.76 

31 .099 

169.49 



A3 

= 188.7 

V 

0.045 

17.48 

4.46 

34. 159 

51.01 

3 

UMSNING 

Al 

= 329.4 

L 

0.136 

41,46 

7.88 

29. 159 

121,41 



A2 

= 328,2 

T 

0. 153 

51 .94 

10, 12 

26.439 

149.39 



A3 

= 264,3 

V 

0.071 

38.31 

12,50 

41 .060 

61 .92 


L : Longitudinal, T : Transverse, V ; Vertical 
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SR. LOCATION 

EP I CENTRAL 

D I RN . 

PEAK 

1 c: v_ L-f I VJ C u 

PEAK 

HI U L L 

PEAK 

arr raoesn 

DURATION 
TR IFUNAC 

NO . 

DISTANCE 


6R0UND 

GROUND 

GROUND 

h BRADY 




ACCL . 

VEL. 

DISPL . 


1 






( 1975 ) 

— 1 

{ km ) • 


< Q ) 

(mm/sec ) 

( mm ) 

(sec ) 


BHATWARI 


lUTTARKASHI 


3 GHANSALI 


RUDRAPRAYAG 


40.2 


54.4 


0.118 

0.117 

0.101 

0.053 

0.052 

0.045 


80.44 

78.21 

95-94 

20.70 

27.06 

17.92 


13.57 

13.37 

25.90 

7.85 

4.01 

3.83 


11.119 
11 . 139 
13.079 

6 . 860 
6.860 
11.56 


Array . 
SPECTRUn 
IIMTENSITY 
FOR 5% 
DAMPING 


0.253 178,73 37.54 13.559 
0.247 297.78 53.23 11.699 
0.294 133.65 23.53 17.879 

0.242 169.56 21.15 15.119 
0.310 194.68 19.85 14.079 
0.196 141.56 22.98 18.560 


128.28 

206.27 
77,09 

207.06 
172. 11 
134.67 

92.52 

91.21 

83.72 

506.97 

579.37 

433.27 


ITEHRI 


SRINAGAR 


BARKOT 


3 KOTESHNAR 


? KARNAPRAYAG 


0 [PUROLA 


0.073 

0.062 


42,15 

92.30 


j 8.17 
19.84 


59.4 


64.6 


65.5 


0.059 

88.41 

23.68 

0.067 

19.45 

5 . 8 b 

0.050 

20.20 

5 . 07 
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f L : Longitudinal, T : Transverse, V : Vertical 
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